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Abstract 



In this thesis, we introduce and evaluate several new procedures for the gener- 
ation and manipulation of entangled states in quantum optical systems. Each is 
evaluated in terms of realistic models of the imperfect apparatus of any real labo- 
ratory implementation. 

First we consider atoms interacting with optical cavities. We describe a method 
of generating maximally entangled states in two separate microwave cavity modes 
by the passage of an atom between them and its subsequent measurement. We 
demonstrate the natural robustness of this approach against a variety of experimen- 
tal imperfections. 

We then introduce a further scheme for the generation of maximally entangled 
states between pairs atoms or ions trapped in separate cavities. The ions are weakly 
driven by an external field and a leaking photon from one cavity is detected, after 
passing through a beam splitter to erase which-path information. In the ideal case, 
this leads to the deterministic generation of high fidelity states, and when more 
realistic conditions are assumed, still allows high fidelity states to be generated in a 
heralded non-deterministic manner. 

Finally, we turn our attention to the manipulation of entangled states of light 
pulses. Entanglement distillation has many applications in quantum information 
science. Gaussian states are the natural states of light pulses produced in laser 
experiments, and standard optical operations preserve the Gaussian form. Entan- 
glement distillation on Gaussian states is impossible with these operations alone. 
We illustrate that the additional element of photon detection enables one to circum- 
vent this restriction, but produces highly non-Gaussian states. We then introduce 
a simple protocol, which we call "Gaussification" , which allows more highly en- 
tangled approximately Gaussian states to be distilled from a non-Gaussian supply. 
Together, these form an entanglement distillation scheme for Gaussian states, whose 
convergence characteristics and realistic operation we analyse. 
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Chapter 1 

Quantum Mechanics and 
Entanglement 

Entanglement is the feature of quantum mechanics which sets it apart most starkly 
from the classical theories that preceded it. It has long been seen as an interesting 
curiosity, but the establishment of quantum information science has given it a new 
role. The goal of quantum information science is to exploit the features of quantum 
mechanics which differ from classical physics, in order to accomplish tasks which 
conventional reasoning, based implicitly on classical physics, would deem difficult or 
even impossible. 

In this thesis, we will introduce several new proposals for the generation and 
manipulation of entanglement. Each proposal is practically motivated and based 
in the physics of real quantum optical systems. We will assess the feasibility of 
each with detailed simulation of the most important imperfections which one would 
expect to encounter in the laboratory. 

In this initial chapter, we will introduce the non-classical properties of entangle- 
ment and discuss its quantification and some of its applications in quantum infor- 
mation science. First, though, we will briefly review the most important quantum 
mechanical concepts, focusing on the features of the theory which have most impor- 
tance for quantum information science. We will assume that the reader has a basic 
knowledge of linear algebra. Good introductions to linear algebra as applicable to 
quantum mechanics can be found, for example, in [Cohen- Tannoudji77, NielsenOO]. 
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1.1 Quantum mechanics 



Our discussion will be restricted to an overview of the most important concepts. 
More detailed introductions can be found in the references just cited, and many 
other quantum mechanics textbooks. 

1.1 Quantum mechanics 

1.1.1 Quantum systems and quantum states 

A quantum system is a number of physical degrees of freedom in a physical object 
or set of objects which is to be described quantum mechanically. For example, this 
can include the position or momentum of a particle, the internal energies of several 
particles or the energy of an electro-magnetic field. 

The quantum state of a system is a mathematical object which represents the 
knowledge we have about the system, and from which all measurable physical quan- 
tities relating to the system can be calculated. 

A special class of quantum states arc called the pure states. Pure states, denoted 
in Dirac's notation by a "ket" , written {ip), represent a ray in a Hilbcrt space Ti, i.e. a 
complex vector space with an inner product. The dimension of is a property of the 
degrees of freedom being described. For example, the state of a spin-half particle 
lives in a two-dimensional Hilbert space, while the Hilbert space for its position 
has infinite dimensions. Due to their simplicity and its importance for quantum 
information science, we will often use a two-dimensional systems for illustration. 
Such a systems is called a quantum bit or qubit, and its basis vectors are labelled 



where |i) are normalised basis states. Thus pure states are sometimes called state 
vectors. The dual space H* also plays an important role. Its elements are called 
"bras" and written {ip \ . The dual of a vector in H represents a linear functional on 
Ti and thus the inner product {ip\(f)) is a complex number while |^) (0| is an operator 



|0) and 

When we express a pure state in terms of a basis in 7i it takes the form 
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on H. 



The pure states are only a very special class of state, whose properties will be 
described below. More general states cannot be described by a simple state vector, 
but require a density matrix [von Neumann27, Landau27] . A density matrix is a 
positive Hermitian operator on Ti. with finite trace, and a vector in the product 
space Ti x Ti.* It is usually denoted p. A basis in 7i x Ti* is formed from the product 
of the set of basis vectors of 7i with the basis of dual space 7i*. Thus a general state 
may be expressed 



The state must be normalisable, i.e. capable of being rescaled so that it has unit 
trace. This is especially important for infinite dimensional states as it means that 
the trace of p must be finite, or in other words p must be trace- class. When working 
with infinite dimensional states it is sometimes useful to formulate "states" which 
are not trace-class. Such states are called improper and cannot exist in nature, but 
can be useful mathematical tools. It is usually convenient to work with normalised 
states. Normalisation is the multiplication of the state with a real factor such that 
the trace is unity. In this thesis we shall assume all states are normalised unless 
otherwise specified. 

1.1.2 Combining and disregarding sub-systems 

Given two quantum systems with state spaces Hi and 0.2, their joint state is de- 
scribed by a state on the tensor product space T^i (8)7^2- The description of the state 
then proceeds as above within this larger space. For example, the joint state of two 
qubits, the first in state |0) and the second in |1) is written |0) This tensor 
product structure has the consequence, however, that there exist states in the joint 
space which did not exist when the systems were described separately. These are 
the entangled states which will be discussed at length below. 

Given a quantum system made up of several sub-systems, it is natural to ask 
how the system should be described when one or more sub-system is out of reach of 
the experimenters. Can a state of the remaining system be described in a consistent 




(1.2) 
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way? The answer is yes, as long as these sub-systems interact at no later time with 
the sub-systems left behind. The procedure to construct this reduced state of the 
remaining sub-systems is simple. One performs a partial trace over the sub-space 
corresponding to the sub-system which is out of reach, 

P=Y1 Pi,j,k,i\i)\j){k\{l\ ^ ■ ^^•^) 

i,j,k,l j i,l 

In fact, the partial trace is the only way of constructing the reduced state which 
satisfies the following physically reasonable demand. Once the lost sub-system is out 
of reach, no operations on it may affect the measurement statistics of the remaining 
sub-systems^. As we shall see below, the partial trace is an important tool for the 
quantification of entanglement in pure states. 

1.1.3 Dynamics 

The evolution of a quantum system is determined by its Hamiltonian H, the energy 
observable for the system, according to the von- Neumann equation 

P^~[H,p], (1.4) 

the density matrix form of the famous Schrodinger equation. Formally integrated, 
this gives, 

p{t) = Uit)pit = 0)Uity , (1.5) 

where U{t) — Texp[—i H{t')dt' /h] and T is the time-ordering operator, which is 

needed since Hamiltonians at different times may not always commute. We see that 

U{t) is a unitary operator since U{t)W{t) = 1. Thus the most general evolution for 

a closed quantum system is unitary evolution. 

This does not apply for open quantum systems, such as an atom interacting with 

surrounding electro-magnetic field modes, where only a sub-system of the unitarily 

evolving system is described. In this case, the evolution of the state is described by 

^It is this condition in fact which means that even though entangled states have apparently 
"non-local" characteristics, no operation is possible which allows instantaneous signalling between 
distant parties, and thus there is no conflict between quantum mechanics and special relativity. 
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a completely positive map. A positive map is a map which preserves the positive na- 
ture of a state and thus transforms it into another physical state. Complete positivity 
demands that 1 (8) p is also positive. An example of an operation which is positive 
and not completely positive is the transpose operation, which takes a state's matrix 
elements {i\p\j) to their diagonally opposite counterparts {j\p\i). Consider the (non- 
normalised) two-qubit state p= |0,0)(0,0| + |0,0)(1, 1| + |1, 1)(0,0| + |1, 1|. If 
we perform the transpose operation on the second qubit alone wc obtain the partially 
transposed density matrix p^^ = |0,0)(0,0| + |0, 1)(1,0| + |1,0)(0, 1| + |1,1)(1,1|. 
This density matrix is no longer positive since it has a negative eigenvalue —1, and 
therefore no longer describes a physical states. Thus positivity alone is not sufficient 
for an operation to be physical. It must be completely positive. 

Any completely positive map may be expressed in Kraus form [KrausTl], 

p^Y.^,pkl (1.6) 

i 

where Ylii-^l-^i — 1- Physically, the Kraus operators Ki of a completely positive 
map on an open system may be obtained if the unitary evolution on the system 
and its environment is known, by performing a partial trace over the state of the 
environment^. Alternatively, given a set of Kraus operators one may always con- 
struct such a unitary operation over the system and an environment, although this 
construction is not unique. 

1.1.4 Measurement of quantum systems 

The traditional way in which measurements on quantum systems are described is 
in terms of observables. Observables are Hermitian operators which correspond to 
physically measurable quantities such as energy, momentum, spin, etc. Any Her- 
mitian operator has a complete set of real eigenvalues corresponding to orthogonal 
eigenspaces. 

A projective measurement of observable O on a system in state p is described 

^This is only possible if the state and environment are initially in a product state. For an 
example of the difficulties of characterising the dynamics of an open system initially entangled 
with its environment see [Stelmachovic01,Hayashi03]. 
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in the following way. We assume that some measuring device exists which returns, 
as a measurement result, one of the eigenvalues of O. Once a result A has been 
registered, the state of the system is transformed by a projector onto the eigenspace 
belonging to eigenvalue A, Pa, i-e. 

p-^P^pP^. (1.7) 

The probability of that measurement outcome occurring is Tr[Pxp]/Tr[p]. In 
quantum information science it has proven beneficial to introduce more generalised 
concepts of measurement as a counterpart to the generalised evolution described by 
equation ()1.6|) . However, these will not play a role in this thesis and will not be 
discussed here. Excellent introductions to generalised measurement exist in [Peres93] 
and [NielsenOO]. 

1.1.5 Von Neumann entropy 

The introduction of measurement outcome probability leads us naturally to a defini- 
tion of entropy for a state. In classical statistical mechanics, the Boltzmann entropy 
of a classical discrete probability distribution pi is defined = —J2iPi^'^&2Pi- 
an analogous way, an entropy for normalised quantum states was introduced by von 
Neumann and now bears his name, 

= -T:[plog2p] = -^^(AaogsA,) , (1.8) 

i 

where Aj are the eigenvalues of p. If the observable which shares the same eigenspaces 
as p is measured, Aj are the probability of outcome i and the von Neumann entropy 
reduces to the Boltzmann entropy of this probability distribution. The von Neumann 
entropy is an extremely important quantity in quantum information science. As well 
as its importance for quantifying entanglement (see below) it replaces Shannon's 
entropy for classical information as the fundamental unit for the quantification of 
information in quantum information theory (see [NielsenOO], chapters 11 and 12 for 
a review). 

The von Neumann entropy sheds more light on the special status of pure states. 
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as these are the only states of zero entropy. Thus only for pure states does a 
projective measurement exist whose outcome can be predicted with certainty. On 
the other hand, density operators that are proportional to the identity have the 
maximal entropy for a system of their dimension. Since mixed states can always be 
written in diagonal form as the convex sum'^ of a set of pure states, 

p=j2Ximm , (1.9) 

i 

a mixed state p is sometimes interpreted as being a probabilistic mixture of this set 
of states, i.e. with probability Aj the system is in the pure state However, it 
is important to recognise that while such a state is indeed indistinguishable from 
a state prepared in this way, this decomposition is not unique. In fact, there are 
infinitely many possible decompositions for any density operator [Schrodinger36] . 
For example, 

I (|0)(0| + |1)(1|) = \ m + |1))((0| + (1|) + (|0) - |1))((0| - (1|)) (1.10) 

Thus one should be wary of assigning special status to any particular pure state 
decomposition of p. 

1.1.6 Complementary observables and the uncertainty principle 

While individual measurement outcomes arc probabilistic, mean or expectation val- 
ues for the operators are straightforwardly calculated in the density matrix formal- 
ism using the following expression, 

(6) = Tr[dp] . (1.11) 

Uncertainty in measurement outcomes can therefore be quantified in terms of the 
variance of O written [A(0)]^, 

[A(0)]2 = (02) _ (0)2 (1.12) 

^By "convex sum" we mean here a sum of normalised operators with positive coefficients which 
sum to unity. 
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and its positive root A(0), often referred to as the "uncertainty" of O. 

The value of the commutator between observables, [Oi,02] = O1O2 — O2O1 



plays an important role in the relationship between the uncertainties of different 
observables for a given state. If two operators commute, an eigenstate of one is an 
eigenstatc of the other, so the system may be prepared in states of zero uncertainty 
for both observables. If a pair of observables do not commute, the converse is true, 
and a system may never be in a state of zero uncertainty for both. Thus there is 
a relationship between the uncertainties of observables which depends upon their 
commutator. This is the general form of the Heisenberg uncertainty relation. 



As we have mentioned above, entanglement is one of the features of quantum me- 
chanics most at odds with a classical world view. In this section we will introduce 
definitions of entanglement, and discuss how it is incompatible with a world- view 
based in classical physics. We will then discuss how entanglement may be quanti- 
fied, before describing some of the important applications of entangled states which 
have been proposed in quantum information science. Since this thesis will focus on 
bi-partite entanglement (between pairs of sub-systems) we shall restrict our discus- 
sion to this simplest case, and not mention more general multi-partite entanglement, 
which has a more complicated structure that is still far from well understood. 

A pure state on Hi (8) 'H2 is entangled when it cannot be written as a product 
state - a tensor product of a pure state on Hi times a pure state on H2- For example, 
the two-qubit state |0)i (8) |0)2 is a product state. Here we have chosen to label the 
states on both qubits and 1 and distinguish between them via the suffixes 1 and 
2. In this thesis we shall often omit the symbol (8) and write |0)i|0)2 = |0)i (g) |0)2. 
Furthermore for compactness we shall often omit the suffixes and differentiate the 
different systems by the order in which they are written, using the concise forms 



A mixed state is entangled when it cannot be written as a separable state, i.e. as 




(1.13) 



1.2 Entanglement 



0,0) = |0)|0) = |0)i|0)2. 
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a convex sum of product states, 



(1.14) 



n,m 



The definitions arc motivated by the concept of Local Operations and Classical 
Communication, LOCC [Werner89]. We assume that the sub-systems we are con- 
sidering are in spatially separated laboratories where general quantum operations 
including measurements may be performed, which are connected to each other via a 
"classical" communications channel such as a telephone connection. The separable 
states are those states which, starting from sub-systems in a product state, can be 
generated by LOCC, the entangled states are those that cannot. Thus entanglement 
can be seen as a kind of non-classical correlation between the sub-systems. 

It is natural to ask whether, given a particular density matrix, it is easy to ascer- 
tain whether the state it describes is separable or not. In fact, algorithms [Doherty02] 
to check this for an arbitrary matrix can take a long time to compute, belonging to 
the class of algorithms known by computer scientists as NP-HARD [Gurvits02] . For- 
tunately, a simple criterion exists which heralds the entanglement present in some 
states, and indeed, for some simple cases, is a necessary and sufficient condition for 
entanglement. This is known as the positivity of partial transpose or PPT criterion, 
and was introduced by Peres [Peres96] and the Horodeckis [Horodecki96] . It is based 
upon the fact that transposition is a positive but not completely positive operation. 
Thus, if we consider a bipartite quantum system in a separable state pscp, and apply 
the transposition operation to only one of the systems (thus partial transpose), the 
partially transposed state has the following form; 



Since transposition is positive, {\ipi) {ipi\)'^ is still a positive operator and corresponds 
to a physical state, and therefore is positive as well. However, this argument 
does not hold when the state is entangled and entangled states may possess a neg- 
ative partial transpose. While possessing a negative partial transpose is a sufficient 
condition for entanglement, it is not in general a necessary condition, and certain 




(1.15) 



n,m 
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states exist which do have positive partial transpose [Horodecki97]. 

Nevertheless, for certain systems of low- dimensionality the breaking of the PPT 
condition is a necessary and sufficient condition for entanglement [Horodecki96] . 
These include pairs of qubits and systems of one qubit and one three-level system 
as well as the bi-partite Gaussian states, an important class of infinite dimensional 
states which will be introduced in chapter ^ 

1.2.1 EPR argument and the Bell inequality 

Now that we have formally defined entanglement, it is time to investigate some of the 
properties which make it incompatible with the expectations of classical physics. We 
will reproduce the argument of Einstein, Podolsky and Rosen (EPR), who, in their 
seminal 1935 paper [Einstein35] argued that quantum mechanics is incompatible 
with some basic intuitions about the physical world. 

EPR argued that physical quantities which can be predicted with certainty must 
be "elements of reality" and for a theory to be "complete" these values must there- 
fore be incorporated into the theoretical description of the system. If quantum 
mechanics is complete in this sense, it follows that non-commuting observables can 
never be simultaneously "elements of reality" as if one property can be predicted 
with certainty, the other cannot. EPR used the following argument to show that 
entanglement in quantum mechanics leads to an apparent contradiction with this 
assertion and that quantum mechanics is an incomplete and thus unsatisfactory 
theory. 

We shall use the version of their argument introduced by Bohm^ [BohmSl] and 
consider the following entangled two-qubit state, 



and assume that the two qubits are located in two spatially separated laboratories; 
one belonging to Alice and one to Bob^. 

*For their example, EPR used infinitely entangled two-mode Gaussian states which are, strictly 
speaking, unphysical non-normalisable states. However, this does not affect the validity of their 
argument. 

^The names Alice and Bob arose in discussions of classical communication theory and have 
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We will introduce a very important family of single qubit observables, the Pauli 
operators. 

ax = |0)(1| + |1)(0| ay = -^|0)(1| +«|1)(0| az = |0)(0| - |1)(1| (1.17) 

Together with the identity operator, the Pauli operators form a basis for two- 
qubit Hermitian operators, i.e., any single qubit operator can be written as a linear 
combination of ax, cry, crz, and 1. The Pauli operators fulfill commutation relations 
of the form [axjCrz] = —Hay and cyclic permutations thereof. 

Imagine that Alice measures the observable cr^ on her qubit, and obtains the 
measurement outcome +1 or -1. Using the measurement postulate outlined in sec- 
tion we see that after the measurement Bob's qubit will be in state |0) or |1) 
depending on Alice's measurement outcome. Thus, now Alice (but not Bob!) knows 
with certainty the outcome of a az measurement on Bob's qubit. In EPR's language 
this value must be an "element of reality" . 

Now imagine that Alice measures ax instead. This leads to Bob's state being 
projected into either state (l/-\/2)(|0) + |1)) or (l/-\/2)(|0) — |1)), the eigenstates 
of ax, depending on the measurement outcome. Now the value of ax on Bob's 
qubit is an "element of reality" . However, as we have seen, if we consider quantum 
mechanics to give a complete description of reality in the EPR sense, the value of 
non-commuting observables can never simultaneously be predicted with certainty. 
Thus there is a contradiction, and EPR conclude that quantum mechanics is not a 
complete description of reality. 

1.2.2 Bell inequalities 

These ideas upon which EPR based their argument were made more precise by John 
Bell. John Bell showed that the implicit assumption behind EPR's argument is that 
physics should posess local realism, that is; 

1. "Realism" - Every measurable quantity should have a definite value regardless 

been adopted whole-heartedly by the quantum information science community, especially in the 
discussion of spatially separated entanglement. Needless to say, they will appear throughout this 
thesis, and can be taken as shorthand meaning two spatially separated experimenters. 
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whether the measurement is made or not. 

2. "Locahty" - These quantities must be features of "local systems" in the rel- 
ativistic sense, thus no property can be shared by two spatially separated 
systems. 

Bell showed that a theory constructed under these two assumptions must fulfill 
certain inequalities, and that quantum mechanics violates these inequalities [Bell64]. 
We will construct such a Bell inequality here. This example is one of a class of such 
inequalities proposed by Clauser, Horne, Shimony and Holt [Clauser69] . 

To derive a Bell inequality one must assume that physics obeys local realism. 
The class of theories for which this is the case are called local hidden variable models. 
Imagine that Alice and Bob each possess a particle upon which they can make two 
different measurements corresponding to two different physical properties which the 
particles possess. Let us label these properties Ai and A2 on Alice's particle and Bi 
and B2 for Bob and assume that each property can have only two values +1 or — 1. 

Now consider the value of the composite quantity, 

AiBi + A2B1 + A2B2 - A1B2 = {Ai + A2)Bi + {A2 - Ai)B2 . (1.18) 

Since Ai and A2 take the values ±1 this expression must equal either ±2Bi or ±2^2. 
Either way the quantity is equal to either +2 or -2. 

Now let us imagine that the particles belong to some ensemble prepared such 
that the probability that {Ai, A2, Bi, B2} = {01,02,^1,^2} is p(oi, 02, &i, ^2) where 
tti and bi take the values +1 and -1. Now we define E[M) as the expectation value 
for the measurement of M, thus 

E{AiBi) + ^(^2^1) + ^(^2^2) - ^(^1^2) 

= + + ^2^2 - ^^2) (1.19) 

= Pa,,a2M,b2iABi + A2Bi + A2B2 - A1B2) < 2 . 

01,02,61,62 

This is an example of a Bell inequality, an inequality that must be satisfied by the 
measurement results of a theory if it fulfils both the realism and locality criteria 
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above. To see that in quantum mechanics Bell inequalities can be violated con- 
sider measurements on the entangled state \iPepr) introduced above. Let Alice's 
measurements Ai and A2 be observables az and ax, and let Bob's be the linear 



the Pauli operators form a basis for single-qubit observables), which one can easily 
confirm have eigenvalues ±1. 



Since {{o'x)a{o'x)b) = {{<7z)a{o-z)b) = —1 and the other expectation values are 
zero, this gives 



so the Bell inequality ()1.19p is violated. 

This shows that quantum mechanics is not consistent with both assumptions 
1) and 2), and if quantum mechanics is believed to be correct, then these basic 
assumptions about physical reality need to be revised. Thus either quantum me- 
chanics must be considered intrinsically non-local or we must accept a theory whose 
measured quantities do not correspond to elements or reahty. 

The Bell inequalities are not merely a philosophical device. They provide a di- 
rect test for this strange and surprising aspect of quantum mechanics since entangled 
states such as \iPepr) can be created in the laboratories and these expectation val- 
ues measured directly. Notable early experiments to attempt this were performed 
by Clauser and Freedman [Freedman72] and Aspect and co-workers [AspectSl, As- 
pect02]. Although there were many loopholes^ in the experiments, which meant 

^Among other loopholes in Aspect and co-workers' experiments, an important one was that, 
unavoidably in the design of the experiment, many of the measurement results were lost. Since 
it is, in principle, possible that the missed results, would have changed these average quantities 
such that no Bell inequality was violated, this is a serious loophole. In fact, a local hidden variable 



combinations Bi = a/1/2(— — ax) and B2 = ^/l/2{az — ax) (remember that 



E{AiBi) + E{A2Bi) + E{A2B2) - ^(^1^2) 




E{AiBi) + E{A2B{) + E{A2B2) - E{A^B2) = 2^2 > 2 



(1.21) 
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that in principle it did not categorically rule out local realism, these and later re- 
sults (see [Zeilinger99] for a historical review) have agreed emphatically with the 
predictions of quantum mechanics. In fact, implementing a loophole-free test of 
Bell inequality violation has proven very technically challenging; however no experi- 
ments have contradicted quantum mechanics, and experiments which claim to close 
all loopholes are being attempted at the moment [Fry02]. 

We must emphasise that entanglement is fundamental to the Bell inequality 
violation since if one considers a quantum mechanics where only separable states 
are allowed, it is clear that the assumptions of local realism can indeed be satisfied 
and no violation will occur. Thus, not only are the EPR debate and the Bell 
inequalities indications that quantum mechanics can lead to results at odds with 
reasoning based on intuition from classical physics, but they emphasise the special 
role of entanglement in this non-classical effect. 

1.2.3 Quantifying entanglement 

It is clear that in some sense state IV'epr) has a high degree of entanglement since 
measurement results on the two particles are completely correlated and yet individ- 
ually completely random, and on the other hand, separable states have no entan- 
glement. Entanglement quantification attempts to provide a number to express the 
degree of entanglement of the state. 

Quantifying bi-partite entanglement in pure states is very simple. If a pure 
bipartite state is not entangled, when we take a partial trace over one sub-system, 
the state of the remaining sub-system will be pure. However, if the state is entangled 
the reduced state is necessarily mixed. It thus makes sense to associate the degree 
of entanglement of the whole bi-partite state with how mixed the reduced state is. 
The measure of mixedness we choose is the von Neumann entropy of the reduced 
state, since, as we shall see below, this gives the entanglement measure a direct 
operational interpretation. The pure state entanglement measure Ey^ can thus be 

model consistent with Aspect and co-workers' measurements has actually been constructed [San- 
tos92]. For a discussion of other loopholes and some healthy scepticism of quantum mechanics 
see [Thompson04]. 
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written, 



£;vN = -TV[p^log2p^] , 



(1.22) 



where p"^ is the reduced state of one of the sub-systems. Since the entropy of a pure 
state is zero, the entropies of both sub-systems of a bi-partite pure state must be 
equal. 

The pure state entanglement has the following important properties, which make 
it suitable as a measure of entanglement [Vedral97]; 

1. It is zero for separable (product) states. 

2. It is invariant under local unitary transformations. 

3. It is not increased on average by local operations with classical communication 

(LOCC) on the two sub-systems. 

A quantity which satisfies these three properties is called an entanglement mono- 
tone [Vedral98, Plenio98a]. Further properties of the pure state entanglement are 
additivity, i.e. the entanglement of the two copies \ip) (g) \ip) is double the entan- 
glement of a single copy {ip), and asymptotic continuity, which means roughly that 
small changes in the state lead to small changes in the entropy. 

Since the reduced density matrix of the sub-systems of \iPepr) is proportional 
to the identity, the entropy of its sub-systems is maximum, and the state can be 
called maximally entangled. The maximally entangled two qubit states have many 
important properties and play a central role in quantum information science. For 
example, since entanglement can not increase under LOCC, these are the only states 
from which any other state can be obtained asymptotically by LOCC alone^. One 
can form a basis from four maximally entangled states, and these states are called 
the Bell states, named in honour of John Bell, 



""The general problem of which pure states may be converted under LOCC to which other pure 
states was solved by Nielsen [Nielsen99]. 




(1-23) 




(|0)|1)±|1)|0). 
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The Bell state basis has many important applications in quantum information sci- 
ence and will feature prominently in the examples below. 

Although E^i^ is a satisfactory entanglement measure for pure states, one must 
consider whether other measures may exist, which may even impose a different order- 
ing of the states. In fact, one can show [Popescu97,Donald02] that the entanglement 
measure for pure states is unique, in the sense that all other additive and continuous 
entanglement monotones are necessarily monotonic functions of -EvN- 

For mixed states, however, the picture becomes much more complicated. Several 
entanglement monotones for mixed state have been proposed (see references within 
[Plenio98a] and more recent proposals such as [Christandl04]). These measures do 
not always give the same ordering of states [VirmaniOO], making the study of mixed 
state entanglement a complicated affair. 

Here we will focus on several important mixed state entanglement monotones, the 
entanglement of distillation, the entanglement cost and the logarithmic negativity. 
The first two of these compare the state to a maximally entangled Bell state in 
different ways. The first asks, given n copies of the state p, what is the maximum 
number m of Bell states which can be made (or distilled®) on average with LOCC? 
The entanglement of distillation is the ratio m/n in the asymptotic limit, where 
n tends to infinity. The entanglement cost [HaydenOl] goes the other way, asking 
how many Bell states would one need (asymptotically) to prepare p with LOCC. 
Both of these measures are equal for pure states and coincide with the pure state 
entanglement -EvN, which gives i?vN a double operational meaning [Bennett96a] . The 
two measures are not, in general, equal for mixed states. 

For an arbitrary state, particularly of high dimension, these quantities can be 
very difficult to calculate. For this reason, the logarithmic negativity has been in- 
troduced [Zyczkowski98,Eisert99a, Vidal02]. This quantity is inspired by the Peres- 
Horodecki PPT separability criterion, described above, and provides a way of quan- 
tifying to how great an extent this criterion is violated. 

The logarithmic negativity is defined 



E^f = log2 IIp^-^IIi , 




Entanglement distillation [Bcnnett96c] will be discussed in more detail in Section 17?^ 



23 



1.3 Entanglement in quantum information science 

where \ \A\\i is the trace-norm of matrix A, defined as Tr[vA^A]. The trace-norm is 
equal to the sum of the absolute value of the eigenvalues, and thus also singular val- 
ues of A. Since singular value decompositions of large matrices may be accomplished 
efficiently JV and Ej^ can readily calculated even for high dimensional systems. The 
logarithmic negativity fulfils the three criteria for an entanglement monotone [Vi- 
dal02, VirmaniOl] . Furthermore, it is an additive quantity and is an upper bound 
to the entanglement of distillation [Vidal02]. It also has an interpretation as the 
entanglement cost of certain states under operations which preserve the positive 
partial transpose of a state [AudenaertOS] . Its tractability makes it very useful in 
the investigation of entanglement in mixed states, and we shall use it throughout 
part mil of this thesis. 

1.3 Entanglement in quantum information science 

As we have seen, measurements on entangled states of spatially separated sub- 
systems allow physicists to test fundamental notions about the nature of the physi- 
cal world. However, quantum information science gives us a new way to think about 
entanglement, namely as a resource to carry out tasks in data processing and com- 
munication which would otherwise be very hard or even impossible. In this section 
we will briefiy outline some of the most important of these proposed applications. 

Here we will focus on qubit entanglement, but most of the applications can 
be generalised to higher dimensions. In particular, we will leave the discussion 
of applications of entangled infinite-dimensional states until section 16.41 In the 
following applications, maximally entangled states play a special role. They are a 
direct resource for quantum teleportation, quantum key distribution 

1.3.1 Quantum key distribution 

The only system of public cryptography which has been proven unbreakable [Shan- 
non49] is the "one-time-pad" [VernamlQ]. Alice and Bob share a secret random 
string called a key. Alice encrypts her message, encoded as a bit-string, by bitwise 
addition modulo two of the key and bit string. She sends this string of data to 
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Bob, and can even broadcast it publicly. Without the key, no information about 
Alice's message can be recovered, but with the key Bob can reconstruct it easily, 
by the same bitwise addition procedure. As long as they only use the key once, the 
one-time-pad cipher is completely unbreakable [Shannon49]. The problem with the 
system is that Alice and Bob need a secure way of generating and sharing the key. 

Quantum key distribution (QKD) solves this problem. By using quantum sys- 
tems, Alice and Bob can generate a key which is completely random, and known to 
no outsiders. The first quantum key distribution was proposed by Bennett and Bras- 
sard [Bennett84] and required no entanglement. Ekert proposed an entanglement- 
based QKD approach [EkertQla], which we shall describe here. 

Alice and Bob share a supply of known maximally entangled states, let us assume 
they are the Bell state | $"'"). Alice and Bob make measurements on the qubits of 
each entangled pair in one of two bases, the eigenbases of <Jz {|0),|1)} and ax 
{(1^)(|0) + |1)), (1\/2)(|0) - |1))}. For each pair they each choose the basis at 
random. They keep a record of each measurement outcome. They then make public 
the measurement bases that they used for each measurement. After they have 
discarded all outcomes where their bases did not match, they compare a subset of the 
remaining data publicly. If the results correlate as expected, they are assured that 
no eavesdropper has intercepted or tampered with their entangled states. They then 
use the remaining data as a secret key. The Ekert protocol was the first protocol to 
explicitly use an entangled state as a resource. Soon afterwards a further application 
employing entangled states was suggested, quantum teleportation. 

1.3.2 Quantum teleportation 

Unlike many other applications in quantum information science, which utilise en- 
tanglement to improve on what is possible with classical physics or with reasoning 
based on classical physics, quantum teleportation achieves a surprising purely quan- 
tum mechanical task. An unknown quantum state is transferred from one system 
to another, via measurements and classical communication alone. Without employ- 
ing entanglement, if Alice has an unknown quantum state and wishes to transfer 
it to her friend Bob at a remote site, the only way to do so is to physically trans- 
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port the quantum system carrying the state. Since measurement of a single copy 
only reveals minimum information about the state and at the same time destroys it, 
no scheme based on measurement and classical communication will allow a faithful 
transfer of the state. 

On the other hand, if Alice and Bob share a maximally entangled state, say 
the Bell state |^'^), this is not the case. They can indeed "teleport" the state via 
measurement and classical communication. Alice makes a two qubit measurement in 
the basis of Bell states in equation ()1.2Hj) . This projects Bob's qubit onto the state 
{(Xxy {<yz)^\4') , where i,j G {0,1} depend upon the measurement outcome. Once 
Alice has broadcast the measurement outcome to Bob, he can correct for these extra 
Pauli transformations and reconstruct {ip)- 

While transferring unknown quantum states does not, at first sight, appear to 
be an application of much utility, in fact teleportation has had immense importance 
for quantum computation which will be described below, both on the conceptual 
level [Gottesman99, Childs04] and practically, e.g. for transferring data between 
remote sites within a quantum computer. 

1.3.3 Other applications of entanglement 

The field of quantum information science has produced many other applications 
of entanglement. Entanglement between spatially separated sites can be used to 
increase classical data transfer rates through so-called "dense- coding" [Bennett92] 
and separated parties can use entangled states to improve the performance in cer- 
tain competitive games [Eisert99b] or distributed computational tasks [Brukner04]. 
Entanglement between systems in a single laboratory also has many applications. 
Entangled states can allow more precise measurements in, for example, frequency 
standards [Wineland92, Huelga97] and lithography [BotoOO], than would otherwise 
be possible. Additionally, entanglement is of fundamental importance in quantum 
computation. 
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1.3.4 Quantum computation 

It has been long known that simulating quantum systems is difficult, as the complex- 
ity of their description grows exponentially with the system size. Feynman was one 
of the ffist to suggest that this complexity could be countered if quantum systems 
were used as an integral part of the simulation devices themselves [Feynman82]. 
An early idea of a quantum computer was proposed by Benioff [BenioffSO], and 
this was developed by Deutsch into the modern idea of a universal quantum com- 
puter [Deutsch85]. David Deutsch postulated that a computer based on the laws of 
quantum mechanics would be so powerful that it could simulate any physical system 
efficiently. The true power of quantum computation was then demonstrated when 
Peter Shor showed that factoring large numbers could be achieved efficiently on a 
quantum computer [Shor94]. The best known classical algorithm increases exponen- 
tially in length with increasing input size. The number of steps in Shor's algorithm 
increases only polynomially giving an enormous "speed-up" over the classical case. 

Entanglement plays an important role in quantum computation. Firstly, Jozsa 
and Linden showed that entangled states must be present in a quantum computer 
at some stage of an algorithm for it to exhibit exponential speed-up compared to 
classical computers [JozsaOS]. Secondly, information must be encoded in entangled 
states if quantum computers are going to be able to operate in the presence of 
noise [Shor96]. Finally, paradigms for quantum computation have been proposed 
based solely on measurements made on certain entangled states, which act as a 
resource for quantum computation [RaussendorfOl, RaussendorfOS]. 

1.4 Contexts and overview of this work 

In this chapter we have seen the power and importance of entangled states, both as 
a fundamental testing ground for quantum mechanics and as a resource for quantum 
information science. Although entanglement occurs almost whenever quantum sys- 
tems interact, the particular quantum states we have identified above as having the 
most practical applications, the maximally entangled states, such as the Bell states 
of equation ()1.23|1 are not observed in nature. There are several reasons for this, 
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the most important is that as these states are pure, they can only exist in systems 
which have no entanglement with any other system, and thus must be isolated to a 
degree which does not occur naturally. 

Controlled generation of entangled states has therefore been an important goal of 
experiments in quantum information science. Close to maximally entangled states 
have been generated in light pulses [Ou88,Kwiat99,Shapiro00] and between pairs of 
atoms and between atoms and cavities [HarocheOS] (an optical or microwave cavity 
is a system of mirrors which confines and concentrates a field). Current schemes 
for entanglement generation in cavity QED^ have several disadvantages which we 
shall describe below in part in particular, requiring a very precise control of 
the evolution of the sub-systems. Other proposals have been made, in particular 
[Bose99] , which require less precise control of the systems, depending instead on the 
detection of photons emitted from the cavities. These proposals, however, are not 
expected to be robust to common experimental errors. In part |n] we will introduce 
schemes for the generation of Bell states in so-called "cavity QED" experiments. 
The schemes we introduce are based on a powerful but simple idea, using weak 
interactions between systems followed by a measurement on one of the systems, and 
this makes them robust against a wide range of experimental imperfections and far 
more suitable for experimental implementation than earlier schemes. 

While the generation of entangled states of light is now very far advanced, ma- 
nipulating these states for quantum information processing has remained a chal- 
lenge. The elements of linear optics, beam splitters and phase-shifters, which we 
will describe in detail in chapter IHl are the most readily available operations for 
the manipulation of light in the laboratory. However, in order to implement, for 
example, quantum computation in light, the linear optical components must be 
augmented by photon counters in a complicated scheme [KnillOl] which is much too 
complicated to be implemented with present technology. Light pulses are quantum 
systems of infinite dimension. They can thus carry so-called continuous variable 
states, which will be introduced in chapter IHl As we shall describe in section 16. 4[ 
many quantum information science applications have been proposed for the contin- 

^Cavity Quantum Electrodynamics (QED) refers to experiments where the hght field of a cavity 
is probed and manipulated through interactions with individual atoms or ions. 
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uous variable domain. However, entanglement distillation for continuous variable 
states, an extremely important application, had so far been unattainable with re- 
alistic operations. In part IIIIl after a detailed introduction of an important class 
of continuous variable states of light, the Gaussian states, we introduce an original 
scheme for the distillation of continuous variable entanglement, a demonstration 
experiment of which could be implemented in a current linear optics laboratory. 

Throughout this thesis, quantum mechanical electro-magnetic radiation plays a 
central role. With this in mind, the second chapter is devoted to a review of the 
quantum mechanical description of light. 
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Chapter 2 

A Quantum Description of Light 



In this chapter, we will show how a quantum theory of electro-magnetic radiation 
may be constructed by "quantisation" of the classical Maxwell equations. Of course, 
as fundamentally different creatures, quantum theories can never be derived from 
their classical counter-parts. Quantisation is rather essentially an ad-hoc process, 
the aim of which is the creation of a quantum mechanical theory, which in certain 
limits matches the classical theory's predictions. The quantum theory of light de- 
scribed here, originally developed by Born, Heisenberg and Jordan [Born26] and 
Dirac [Dirac27], has been verified by a large number of experiments, many of which, 
especially in recent years, have explicitly tested the theory's non-classical predic- 
tions. The treatment presented here is slightly more general than that found in 
textbooks such as [Knight83, Mandel95, Meystre99, LoudonOO] as we shall explain 
below. 
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2.1 Classical electro- magnetic waves 

Classically, the electro-magnetic field in a vacuum is completely described by the 
in-vacuo Maxwell equations [Maxwell73] , 

VxE = -— (2.1a) 
dE 

V X B = /ioeo^ (2.1b) 

V-E = (2.1c) 
V-B = (2. Id) 

where E and B are electric field and magnetic induction and /io and eo are the 
permeability and permittivity of free space. It is convenient to introduce the vector 
potential A which is defined as 

B = V X A (2.2) 
and choose a convenient gauge in which to work, the Coulomb gauge, where 

V-A = (2.3) 

and 

E^-§. (2.4) 

Since there are no external sources, the scalar potential can be set to zero without 
loss of generality. Substituting these equations into the second Maxwell equation 
()2.1bj) and employing the identity VxVxA = V(V-A) — V^A, results in the 
following wave equation for A, 

V^A=^-^ (2.5) 



where c = A/l/(/ioeo) is the speed of light. This wave equation can be solved by the 
method of separation of variables which, in the spatial case, leads to the Helmholtz 
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equation, 



V\(r) 




c- 



(2.6) 



where, for later convenience we choose the separation constant — cj^/c^ and where 
k is a yet unspecified labeL The solutions to equation ()2.6|) . Uk{r) are called mode 
functions and depend on the boundary conditions. Usually, when the quantisation 
of the electro-magnetic field is presented one usually considers only the simplest case 
of mode-functions that are one-dimensional and sinusoidal. We take a more general 
approach in order to show that the same quantisation method is valid for any set of 
mode- functions which satisfy equation ()2.(ij) . and shall thus leave the mode-functions 
unspecified. The label k can take continuous or discrete values, but here we consider 
the discrete case which applies, for example, to the field in a cavity. A discussion of 
the continuous case may be found in, e.g. [Mandel95, page 512]. 

In classical optics, the mode functions represent resonant modes, and in the 
quantum case provide the basis for the quantisation of the field. We choose an 
orthonormal set of mode functions^ and can thus expand A as the sum 



where we have guaranteed that A(r,t) remains real, while allowing complex mode 
functions which are often convenient for calculations. Substituting this mode ex- 
panded solution back into equation ()2.5|1 produces the following equations of motion 
for the coefficients Ak{t); 



This equation is similar to the equation of motion for a simple harmonic oscillator 
and as we will shortly see, the quantisation of the electro-magnetic field is based 
upon this correspondence. Choosing positive frequency solutions of equation ()2.8j) ^. 
and employing the unit polarisation vector e^, Afc(t) = Ok exp[— icutjek we write the 

^This will always exist as is a Hermitian operator on a given region of M'^ under the condition 
that functions and their first derivatives vanish on its boundary, as in this case. 

^This does not restrict generality, because the mode functions always exist in complex conjugate 
pairs the negative solutions will appear in association with the complex conjugate of Uk. 





k 



Afc(t) 



^ dt^ 
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general solution to equation ()2.5j) 

A(r, t)=tJ2 («ke-^"'^*Mk(r)ek - a;;e+^-''*<(r)e^;) . (2.9) 

k 

Thus, from equations ()2.2|1 and ()2.4j) . we obtain solutions for the classical fields E 
E(r, t) = J2^^ (ake-*'^''V(r)ek + 4e+*-''*<(r)e*) , (2.10) 

k 

and B 

B(r,t) = (ake-*-''*(V x (Mk(r)ek)) - aj;e+^-''*<(r) (V x {K{r)el))) . (2.11) 



The classical Hamiltonian of the electro-magnetic field is 



where the integration is over the whole space as defined by the boundary conditions. 
We substitute the mode expansions for E and B (equations ()2.10|) and 1)2.111) ) into 
this Hamiltonian. After some manipulation^, this integral can be solved by using 
the orthonormality of the mode functions which generates the following expression, 

H = 2eoJ2^l\^k\^ ■ (2-13) 

k 

To show that this is equivalent to the Hamiltonian for a simple harmonic oscil- 
lator associated with each mode, let us make the linear transformation of variables. 



a^^J±(x^ + z^) (2.14) 



4eo V (jj 



a* = W-^ fxk-2— ) (2.15) 



4eo V uj 



•^The integral over B is brought into a suitable form in several steps, first employing integration 
by parts and the identity (A x V) • (V x B) = A • (V x V x B). Then the identity introduced 
above allows one to use the zero divergence of Wk(r)ek and the Helmholtz function to remove all 
differential operators from within the integral. 
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in terms of which the classical field Hamiltonian is 

^=^E(^^k + ^k4) ■ (2.16) 

k 

Thus each mode of the field can be associated with a simple harmonic oscillator 
with position and momentum p^. The quantisation of the field then proceeds by 
replacing these classical variables with non-commuting quantum operators, whose 
commutators are specified by the Poisson brackets of the classical theory. 



2.2 The quantum harmonic oscillator 

The Hamiltonian for the quantum harmonic oscillator with unit mass and spring 
constant o;^ is 

H^^-[p'+u;'x'], (2.17) 

where x and p represent position and momentum respectively and obey the commu- 
tation relation [x, p] = ih. The Hamiltonian can be diagonalised by introducing the 

ladder operators; 

a = y ^ (c^x + ip) (2.18) 



a^ = y^(^£-^p) (2.19) 

with commutation relation [a, a^] — 1. Written in terms of the ladder operators, the 
Hamiltonian takes on a simple form 

H^huj (^a^a + . (2.20) 

The eigenstates of this Hamiltonian are called Fock states. Their properties can be 
derived as follows. Let |0) be an energy eigenstate with energy 

H\ct>) = fiw {^^a + 10) = E^\ci>). (2.21) 
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By applying the operator a to both sides of this equation, and using the com- 
mutation relation, we find 

Ha\(l)) = (E^ - huj)a\(t)) . (2.22) 

Thus, the state a|0) is itself an (as yet unnormalised) energy eigenstate with a lower 
energy — fiui. This is why a is often called a lowering operator and repeated 
application of a will lead to lower and lower energy eigenstates. 

That this lowering cannot go on indefinitely, i.e. there is a ground state of lowest 
energy, can be seen by remembering that any state must have a non-negative norm, 
including the state a|0). Therefore (0|a'^a|0) > 0, which bounds the expected energy 
from below, {n\H\n) > hw/2. Thus there exists (at least one) state of lowest energy, 
which we will label |0). Since action of a cannot lead to another physical state, 
this state must fulfil a|0) = 0. This allows us to calculate its energy, and we find 
^|0) = huj/2. 

Just as operator a is a lowering operator, its conjugate a) is a raising operator 
which maps one energy eigenstate to another with energy raised by TvjJ. Thus, 
applying d) n times to the ground state is {l/2)fiuj generates an energy eigenstate 
of energy hw{n-\- {1/2)). We label this eigenstate 

H\n) = huj + \n) (2.23) 

and thus 

a)a\n) = n\n) . (2.24) 

We thus name a^a the number operator and label it n. So far we have not 
considered the normalisation of |n), since the eigenvalue equations arc invariant 
when the state is multiplied by a scalar. We write the state \n) = a„a^^|n — 1), 
where a^n is a normalisation constant which will be determined inductively. If we 
demand that \n) has unit norm, 

|o;n|^(?T' — l|aa^|n — 1) = |Q;n|^(n — 1|(7t, -|- l)|n — 1) = |Q;„|^n = 1 . (2.25) 
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Thus a„ = 1/ ^/n and by induction 

\n) = \l^{dr\0) . (2.26) 
V nl 

It also follows that 

a\n) = y/n\n ~ 1) (2.27) 

and 

a^n) = Vn+ l\n + l). (2.28) 

To show that the vectors \n) form the basis of an infinite-dimensional Hilbert 
space, one still needs to show that each eigenstate is non-degenerate. This can 
be proven under the assumption that only analytic functions of d and act on 
the space, since it can then be shown that the only operators on the space which 
commute with n are functions of n and therefore the eigenstates of n must be non- 
degenerate (see [Messiah72] volume I, page 436 and Problem XII. 1, page 460). 



2.3 The quantum light field 



The quantised version of the electro- magnetic field Hamiltonian, equation ()2.16j) . is 
thus 



k 



2 -2 \ 

2 ' '^k^kJ 

^ (2.29) 



Since the contribution of the sum of ground state energies to this Hamiltonian is 
infinite, we shall adopt here the conventional practise of omitting it and redefining H, 

H = Y^ hukdldk . (2.30) 

k 

An infinite dimensional Hilbert space is associated with each mode, and the state 
of the field as a whole inhabits an infinite tensor product of these spaces. By com- 
paring the quantum Hamiltonian in equation ()2.30|) with the classical Hamiltonian 
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in equation ()2.13j) . we see that the quantum case can be obtained from the classical 
by replacing coefficients with the operators A/^/(2eou;k)ak. Thus by performing 
the same substitution in equation ()2.1U|) we can write a quantum electro-magnetic 
field operator E, 

E = ^ £k (akMk(r)ek + aj^M* (r)e*) (2.31) 

k 

where = A/(^tUk)/(2eo). 

Note that, unlike its classical counterpart, this operator is time- independent. 
This is because we are working in the Schrodinger representation where this time 
dependence is carried by the state. 

2.4 Quantum states of a light mode 

States of a single light mode of a quantum harmonic oscillator are density matrices 
on an infinite-dimensional Hilbert space 

p = Y^apJp){q\ . (2.32) 

p,q 

Since full description of the state requires specification of an infinite number of 
coefficients, it is often convenient to express states in terms of functions in phase 
space. Rather than working directly in the phase space defined by the position 
and momentum variables x and p, we shall use the quantum optics convention of 
defining dimensionless and rescaled operators X = ^Jojjfix = ^/l/2{d + a^) and 
P = a/ / {huj)p = —i^/l/2{a — a'^). These operators are known in quantum optics as 
quadrature operators'^ and obey the commutation relation [X, P] = i. 

A variety of phase-space functions are used in the literature, and each may be 
obtained from the other by the appropriate transform. However, in this thesis we 
will restrict our discussion to the Weyl function or characteristic function of the 
state. 

■^Note that the precise definition and normalisation of the quadrature operators is a matter of 
convention and, unfortunately, there is httle uniformity in the hterature. For a comprehensive 
discussion of all the different degrees of freedom one has in these conventions see [KriigerOl], 
section A-2. 
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First we define the Weyl operator W{X, P) as follows, 

W{X, P) = e^(^^-^^) (2.33) 

where X and P are real valued variables which represent a point in the rescaled 
position and momentum phase space. Apart from a change in the argument, the 
Weyl operator is the same as the phase space displacement operator introduced by 
Glauber [Glauber63], Da, where a = ^JXpliX + iP\ and we shall employ several 
of the properties derived by Glauber in part lllll 

The characteristic function of the state p, xp{?^, P) is then defined 

Xp(X,P) = Tr(pW^(X,P)) . (2.34) 

This function is uniquely defined by and is in one-to-one correspondence to the 
state p and therefore provides an alternative means of representing the state. The 
characteristic function plays an important role in the description of the Gaussian 
States, an important and widely occurring class of state in quantum optical systems, 
which will be discussed in more detail in part lllll of this thesis. 



2.5 Coherent states 

An important class of single-mode states are the coherent states. These states are 
parametrised by a single complex number a as follows. 



a 

They have many important properties, but in particular they behave in many 
ways most closely to classical light. For example, the expectation value of the 
operator a, (a|a|a;) = a. Since under free evolution, these states evolve in time 
as \ipit)) = |e^*'^*a;), it is easy to verify that the expectation value of the electro- 
magnetic field for these states fulfils the equation of motion for the classical field 
mode in equation ()2.10|) . Additionally, one can transform the states and operators 
of the mode using the displacement operator D^a- In this new picture, the state is 
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transformed to the vacuum state and field operators are now tfie sum of a complex 
number a and a bosonic operator. Since this operator acts on the vacuum, its 
effect can often be neglected, especially when \a\ is large, and the mode can be 
represented to good approximation by the "classical" complex amplitude alone. In 
such a scenario, the mode is often described as "classical" . 

2.6 The quantum state of laser light 

A detailed description of the generation of laser hght is beyond the scope of this 
introduction, and many textbooks, for example [Sargent 74] or [Mandel95] provide a 
clear and detailed treatment. Here we will restrict ourselves to commenting on how 
laser may be described quantum mechanically. In the most general of terms, lasers 
work via the following principle. Laser light is generated by stimulated emission 
of light from material in a cavity resonator (see the next chapter for a detailed 
discussion of optical cavities) which is continuously re-excited so that the stimulated 
emission carries on. The resonator typically has one semi-silvered mirror so that the 
cavity radiation leaks out to form a beam of light. 

Under ideal conditions the state of the cavity mode is a Poissonian mixture of 
states of Fock states 

p = e-^^^— |n)H (2.36) 

n 

where A is a real parameter. However, this may be re-written using the following 
identity, 

e-^' J2 = ^ / d(l)\Xe'^){Xe"^\ • (2.37) 

n 

where \Xe'^'^) is a coherent state. What this means is that the state of the cavity is 
indistinguishable from a coherent state of known amplitude but completely unknown 
phase^. 

''This is not, however, to ascribe either the latter or former interpretation of the state. Both 
interpretations are equally valid. 
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To a first approximation, the state of the hght-beam which emerges from the 
cavity, may then be written [Blow90,van Enk02] 



where the hght beam is treated as a sequence of wave-packets, each with a travel- 
hng mode function u{x, y)S{z — ct, w) where we assume that the beam is travelhng 
along the z-axis. The function u{x,y) describes the transverse profile of the beam, 
while S{z,w) = 1 when w/2 < z < w/2 describes a sawtooth function of width 
w [van Enk02]. The width w of the wave-packets can take any value but must be 
much greater than both the (mean) wavelength of the light and the size of the laser 
cavity. 

The correlation in the phase between different wave-packets here is called phase 
coherence. In any real laser, noise in the lasing process means that the phase coher- 
ence does not last throughout the entire beam. Rather, the phase becomes gradually 
decorrelated such that beyond the coherence length separated wave-packets become 
completely decorrelated. 

In part IIIII of this thesis, we will consider the manipulation of the quantum 
state of light pulses. Laser pulses of femtosecond duration are generated in the 
laboratory from a mode-locked broadband laser such as a Titanium Sapphire laser 
(see [Krausz04] for the latest developments in the generation of ultra-fast pulses). 
The laser pulses can be treated as individual wave-packets with an independent 
state space and, ignoring phase decorrelation effects, the state of the pulses may be 
written; 



where we now quite naturally divide the beam to wave-packets consisting of individ- 
ual pulses. Since current pulsed laser sources can have a coherence length of up to 
many 10s of metres [Silberhorn04] and pulses have a typical separation on the order 
of a metre, pulse trains of many pulses can still retain excellent phase correlation, 
as was demonstrated for neighbouring pulses in a recent experiment [Glockl04]. 
In quantum optical experiments, one typically splits the beam at its source with a 




(2.38) 




(2.39) 



40 



2. 6 The quantum state of laser light 




aout 



bout 



Figure 2.1: A schematic diagram of the action of a beam spUtter on a pair of input 
modes, a and b. 

beam-splitter into an "object beam" , that is to be probed and manipulated further, 
and a phase reference beam. 

A beam splitter is a semi-silvered mirror which mixes two incident light beams 
together. Under ideal circumstances, when the two modes are perfectly matched 
in position, timing, frequency and mode envelope, it has a simple action, namely 
that the mode operators, and are linearly transformed to (a''")' = Ta^ + Rb\ 
(fet)' = -i?*at + Tb\ where {T,R} G C and \T\^ + \R\'^ = 1. The magnitude and 
phase of the transmittivity T and reflectivity R are properties of the beam splitter. 
For a more detailed discussion see section lfi.2.1[ When a coherent state lAe**^) is 
mixed with a vacuum mode on a beam splitter the output modes are a tensor 
product of coherent states iTAe*"^) ® l-RAe**^) with a well-defined relative phase. This 
means that the state of a split laser pulse may be written, 



This state is indistinguishable from two phase correlated pulses of unknown 
phase. It is impossible to measure the phase of a coherent state absolutely. It 
has to be measured relative to some reference beam. When one uses the reference 
pulse for all phase measurements, measurements will be consistent with the object 
pulse having been in a coherent state. 




(2.40) 
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In fact, the correspondence goes further. Imagine that some quantum operation 
with Kraus operators Ki is applied to the object pulse 

J2Kip,pmpuis.Kj = i- /" d<f)[J2K,\TXe'^){TXe''^\Kj] ® \RXe"^){RXe'^\ . (2.41) 

When the phase of the object pulse is measured relative to the reference beam, 
the measurement outcome will be indistinguishable from measurement on the state 
[Y.,Ki\TXe'^){TXe''^\Kl]. Since the Kraus operator formalism includes as a special 
case projective measurements, this means that as long as all phase measurements 
on the object beam are performed in relation to the reference beam from the same 
laser source, the state of the beam behaves as if it were a coherent state [M0lmer97]. 

For this reason, it is standard practice to treat the light generated from a laser 
source as a coherent state, as this usually simplifies both the mathematical analysis 
and the conceptual interpretation of experiments. This is the convention we shall 
adopt here, and we shall thus write, for example, |Ae*^) when we mean 

J d(l)\Xe'^'^+'^^){Xe'^^+^^\ (g) |Ae''^)ref(Ae''^| . (2.42) 

The presence of the reference beam will always be implicit. We will introduce 
phase decorrelation as a source of noise when we consider imperfect experimental 
conditions in chapter [8. 4[ 

This is the end of part |lj the general introductory section of this thesis. In 
partHH we focus on entanglement generation with atoms and cavities. In part lllH we 
shall consider manipulating the entanglement properties of light pulses, in particular 
introducing a scheme for entanglement distillation. 
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Part II 

Robust Entanglement Generation 
in Atoms and Cavities 
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Chapter 3 

Cavity QED and Atom-Field 
Interaction 



In this part of the thesis, we shall present two novel approaches for generating 
entanglement in the domain of cavity quantum electrodynamics or cavity QED. 
Cavity QED is the name given to experiments which probe the quantum nature of 
the electro-magnetic field trapped in a cavity, which usually consists of two or more 
mirrors. These mirrors can be of an extremely high quality, so quantum states of 
the resonant cavity modes with lifetimes up to a millisecond [RaimondOl] can be 
observed. Additionally, single atoms [Ye99] and ions [KellerOS] can be trapped inside 
the cavities, such that their interactions with the cavity modes can be manipulated 
and observed. 

In this chapter we shall review the basic theory of cavity QED, including descrip- 
tion of the interaction between an atom and a mode of radiation, especially a cavity 
mode. The methods and results presented here were collected from a number of 
sources, especially [Mandel95], and also including [KnightSO], [Cohen- Tannoudj 189] 
and [LoudonOO]. 

3.1 Cavity modes 

An optical cavity is a system of mirrors confining an electro-magnetic field. The 
mode functions of the cavity are solutions of the Helmholtz equation ()2.6|) under 
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boundary conditions defined by the geometry of tlie cavity mirrors. Many exper- 
iments in cavity QED employ pairs of spherical mirrors, in what is known as a 
Fabry- Perot configuration, and it is this setting on which this discussion will centre. 

Approximate mode-functions for a Fabry-Perot cavity were first calculated nu- 
merically by Fox and Li in [Fox61] in 1961, and approximate analytic solutions 
soon followed, a comprehensive review of which was presented by Kogelnik and Li 
in 1966 [Kogelnik66] . For circular spherical mirrors in a Fabry-Perot configuration, 
the mode-functions are Laguerre- Gauss functions (products of Laguerre polynomials 
and Gaussians). The exact form of these mode functions are given in [Kogelnik66] . 
The modes which are most often coupled to atoms in cavity QED experiments are 
known as TEM„oo modes, the transverse electric and magnetic modes which have 
a Gaussian spatial profile. For two spherical mirrors sharing a central axis and 
with radius of curvature R whose centres are separated by distance L, the TEM„oo 
mode-function is 



where cylindrical coordinates are employed. The z-axis lies along the line connecting 
the centres of the two mirrors and the origin is the centre of the cavity, as illustrated 
in figure 13.11 

The function $(r, z) is 



cos[<l>(r, z)]e '^^w^ 



(3.1) 




(3.2) 



where the wave-number k is 




(3.3) 



The function w{z) is the cavity beam radius 




(3.4) 
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Figure 3.1: The layout of a Fabry-Perot cavity with spherical mirrors of radius of cur- 
vature R separated by distance L. 



The frequency of the mode u is given hj u = ck/{2TT) and we can see from 
equation that the frequency difference between TEM^oo modes will be c/(2L). 
For a microwave cavity where L is on the order of mm, this frequency difference is on 
the order of GHz or THz^. In the literature (see for example [Thompson92,Hood98]), 
^The frequency difference between non-Gaussian modes for a given principle mode number n is 



much smaller than this, actually of the order (c/tt) arctan (^L/{y/L(2R — L))j (see [Kogelnik66] ) 
and there is degeneracy in the non-Gaussian modes. However, this frequency difference is still on 
the order of GHz. That the modes are so widely spaced in frequency will allow us to simplify 
significantly the description of an interaction with an atom which is resonant with a particular 
mode - see section rOl 



with 




\^L{2R-L) 
2^r 



(3.5) 
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the following approximation mode function is often employed 



UnOO 



e ™(' cos{kz). 



(3.6) 



Since this eases calculation considerably and this is a good approximation when the 
atom remains close to the plane at z = 0, we shall use this in our calculations in 
chapter |3 From basic symmetry arguments, it is clear that the mode polarisation 
vectors must he perpendicular to z-axis. 

3.2 The interaction between an atom and a light 
mode 

A full treatment of the interaction between an atom and a mode of the electro- 
magnetic field (see, for example [Mandel95], chapter 15 or [Cohen- Tannoudj 189]) is 
rather complex and beyond the scope of this thesis. For the parameter ranges of 
cavity QED, one can employ the dipole approximation which simplifies the descrip- 
tion of the interaction significantly. This approximation, which treats the atom as 
a point-like dipole, is valid when the wavelength of the light is much larger than the 
atomic dimension and when the light is not of high intensity. 

The classical Hamiltonian for the interaction between a dipole with dipole mo- 
ment D and a homogeneous electric field is 



This Hamiltonian can be cast in quantum mechanical form by replacing D and 
E with the corresponding quantum mechanical operator. The quantum electric field 
operator E has been derived above. We simplify the analysis further by adopting 
a two-level approximation for the atom, where we assume that the atom has only 
two internal states, 1^4) and \B) with energy difference Tiujab- A justification for this 
approximation is that, as we shall see explicitly below, energy exchange between 
the atom and cavity only occurs when the atomic transition frequency is close to 
resonance with the mode frequency. Interaction with non-resonant modes leads to 



Hi = D.E(t) . 



(3.7) 
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energy shifts, which can be accounted for in the way we define the energies of levels 
\A) and \B). In cavity QED experiments, the atomic species are chosen such that 
only the desired transitions fulfill the resonance requirements and thus the two-level 
approximation is valid^. 

In this two-level picture , the most general form that the dipole operator D can 
take is 

D = Baa\A){A\ + Bbb\B){B\ + Bab\A){B\ + Bba\B){A\ (3.8) 

where Djj = (2|D|j). 

The dipole moment operator is proportional to the relative position vector be- 
tween the nucleus and the interacting electron D = — ef and thus has odd parity 
(it changes sign when r is replaced by — r) . Since diagonal components of the inter- 
action operator merely shift the energies of the system's basis states we can neglect 
the diagonal components and set T)aa and D^s to zero. Furthermore, since D is a 
Hermitian operator, D^^ = D^^. For brevity of notation, let us label D^^ = d. 

The interaction Hamiltonian can then be written 

Hi = ^£k((d- eK(r)ak|5)(A| + (d* • e*)ul{ral\A){B\ + 

(d* ■ e)uUr)a^\A){B\ + (d ■ e*)n* (r)aj,|5)(A| 

The analysis of this interaction Hamiltonian can be simplified greatly by the 
use of the Rotating Wave Approximation (RWA). This is another approximation 
which works very well for interactions with non-intense fields. Under the RWA, the 
latter two terms of equation ()3.9j) are omitted. These are the terms which contain 
products of terms which both create or both destroy quanta of energy. It can be 
shown that the higher order effects to which these terms contribute are only non- 
negligible in very intense laser fields (for an example of this, see [BrowneOO]). The 
term "Rotating Wave Approximation" is due to the fact that, when an interaction 

^Sometimes more complicated situations arise where interactions between three, four or more 
levels must be taken into account. As long as each transition is resonant with just one mode then 
this can be treated as a sum of individual two-level interaction Hamiltonians, as is the case in 
chapter |S| 



48 



3.3 Interaction with a single cavity mode 



picture (see section 13 .(ij) is adopted, the neglected terms rotate at a much higher 
frequency than the remaining terms, which means that they have a much smaller 
effect on the system's evolution. 

The interaction Hamiltonian under the RWA takes the concise form 



where gk{r) = {d-e)uk{r)/h contains all spatial dependence of the coupling strength. 
This is often called the Jaynes-Cummings Hamiltonian, after the authors [Jaynes63] 
who first introduced it (see [Shore93] for a review). 

3.3 Interaction with a single cavity mode 

To better understand how the atom will interact with the whole electro-magnetic 
field, let us first consider the interaction with just a single mode. We shall see 
that the behaviour of this interaction depends greatly on whether the energy of the 
transition between the atomic levels is at resonance with the mode frequency. 

As we are only considering a single mode, let us drop mode labels and write the 
full Hamiltonian of the atom-mode system, under the approximations introduced 
above. 



H = hua^a + -^{\B){B\ - \A){A\) + h{ga\B){A\ + g*a^A){B\) , (3.11) 

where u is the mode frequency and hcoAB is the transition energy between levels \A) 
and l-B). 



When dealing with resonance phenomena it is useful to define the detuning, here 
S = uj — uoabi between the transition frequency and mode frequency. Hence, we 
rewrite equation ()3.11|) 



Hi = hY, [g,.{v)a^\B){A\+ gl{v)ai\A){B\) 



(3.10) 



k 



H = huja)a + 



2 



{\B){B\-\A){A\) + h{ga\B){A\+g*a 



\A){B\) . (3.12) 
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Due to the RWA, this Hamiltonian only couples states in pairs \B,n) and | A, n + 
1) , and the Hamiltonian can be rewritten as a direct sum of operators Hn which act 
only in these sub-spaces, 



H = J2Hn, (3.13) 



n=0 

where Hn can be written in matrix form as follows, 



( \A,n+l) \B,n) ) 

h 



uj{n + l) + 5/2 g*Vri + l 
gVn + 1 uj{n + ^)-S/2 




(3.14) 



Diagonalising this matrix for all n will give us the eigenbasis of the full atom- 
cavity Hamiltonian. Let us label these eigenstates, which are often called dressed 
states |a+,n) and \a_,n) [Jaynes63, Cohen- Tannoudji69]. Their respective energies 



arc ?i{ijj{n + |) ± h{ytn/2) where Vtn = a/^^ + 4|5fP(n + 1) is known as the Rahi 
frequency. 

The states themselves can be written 

|Q;_,n) =sin(e)|A,n+l) + e'^''s(^^cos(e)|S,n) (3.15) 



|a+,n) = -cos{d)\A,n+l) + e''''^^^^sin{9)\B,n) (3.16) 



where 



tan(^) = . (3.17) 



3.3.1 Far-detuned behaviour 

Let us first consider the case of far detuning, i.e. when \5\ 3> g. This means that 

tan(^) is approximately zero and therefore |a_,n) \B^n) and \a^,n) ^ |A, n + 1). 
This can be interpreted as follows; the product basis states 1-8,^) and + 1) are 
approximate eigenstates of the complete Hamiltonian, therefore a system prepared 
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in one of these states will remain there. In other words, at far off-resonance, no 
energy exchange occurs between the atom and the cavity mode. However, states 
\a±,n) do not have the same energy as j^l, 1) and \B,n). Thus, the observed 
energy levels of the atom appear to have shifted from their respective natural values 
by =F^(^ ~ ^n)/2 ~ ^Tig^{n + This effect is known as the light shift. 

3.3.2 Near resonance 

Now let us consider the opposite extreme, when \5\ <^ (?. In this limit tan(^) is 
approximately unity and the dressed states are approximate equal superpositions of 
the states \A, n + 1) and \B, n). 

Consider a system prepared in the product state When we expand this 

state in terms of the dressed states its time evolution is easy to calculate, 

\'4^B,n{t)) = e-^'^'-s(«)e-^'^("+^)* (^sin(^)e-*^*|a+,„) + cos(^)e+^^*|Q;_,„)) 



-i^sm \—t ] \A,n+ 1) 



+ 



cos 



s . 

t I + I— sm 



\B,n) . 



(3.18) 



Likewise, the evolution of a system prepared in state 1^4, n + 1) is 



cosl-tU.-sm -t 



\A,n+l) 



.2^ • /^^"+ no \ 



(3.19) 



In the hmit that \6\/g tends to zero, these become 



-iu>{n+^)t 



t ) 1) + COS ( %t ) \B,n) 



(3.20) 
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\i^A,n+l{t)) = e 



■iaj(ri+i)t 




,« arg g 



sm 




jt) \B,n) 



) 



(3.21) 



The atomic population oscillates between ground and excited states with fre- 
quency Qn- This phenomenon is known as Rabi oscillation, which is why ^7„ is 
called the Rabi frequency. 



3.4 Mult i- mode case - interaction with a resonant 



In a standard atom cavity experiment, an atom is prepared in a state with a transi- 
tion close to resonance with a particular cavity mode. Typically, the Rabi frequency 
will be on the order of kHz, while the difference between mode frequencies will be 
on the order of GHz or THz. This means that the far-detuned limit is valid for all 
modes except those near to resonance. 

Thus, while a full treatment of the evolution of the system would require the 
full multi-mode Hamiltonian, it is usually a very good approximation to neglect the 
interaction with the non-resonant modes, since they only contribute level shifts to 
the system's evolution. The single-mode Hamiltonian, with the off-resonant energy 
shifts included in the state energies, describes the evolution of the system very well 
and is the starting point for many research papers on cavity QED. 

3.5 Interaction with a classical field 

When a laser mode interacting with an atom with a Jaynes-Cummings Hamiltonian 
is in a coherent state with large amplitude solely quantum effects, such as 
entanglement between the mode and the atom become negligible, and the state of 
the mode can be treated "classically" to a good approximation. As discussed in 
section 12.51 this involves assuming that the field is in a coherent state |a) with 
I a I ^ 1. One then replaces the mode operators in the Hamiltonian by their time- 



mode 
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dependent expectation values, neglecting any effects on the time evolution from the 
interaction, e.g., operator a is replaced by the expectation value {ae~'^^*\a\ae~'^^*) — 
Q,g-«wt The free field energy in the Hamiltonian is now just a number, and can be 
neglected, as it just imparts a non-detectable global phase to the systems evolution. 
The semi- classical (as the atom is still treated quantum mechanically) atom-mode 
interaction Hamiltonian is 

i^int: S.C = gae-'^'\B){A\ + h.c. . (3.22) 

Typically, the phase of ga is unimportant, so g is assumed to be real, and the 
classical Rabi frequency = 2giQ; is introduced, to give the standard semi-classical 
interaction Hamiltonian under the rotating wave approximation, 

^int: S.C = ^e-'^'\B){A\ + h.c. . (3.23) 

3.6 Simplifying the Hamiltonian: the interaction 
picture 

In quantum optics, the most interesting part of a systems evolution is often the part 
caused by its interactions. The evolution of the system due to the free energies of 
the sub-systems is usually of little interest, and only complicates calculations. For 
this reason, one often employs the interaction picture. 

Let us define Hi^t as that part of the Hamiltonian due to the interaction one is 
interested and Hq — H — H-^^^ as the non-interaction part. Without the interaction, 
the system would evolve via the unitary U{t) — exp[— iifo^/^l- interaction 
picture, one works with the state without this free evolution |'0i.p) = U\t)\%l)) and 
operators are correspondingly obtained as O/.p = U{t)Oi_pU\t). 

Since the Hamiltonian determines the time evolution of the system and U{t) is 
time-dependent, the interaction picture Hamiltonian must be obtained in a different 
way. By inspecting the Schrodinger equation, one sees that the interaction picture 
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Figure 3.2: In a real cavity, the mirrors will not be perfect, photons will leak through 
the mirrors with rate 2k. Also, the cavity will not be totally sealed from the continuum 
of vacuum modes of the surrounding environment. An atom in the cavity will undergo 
spontaneous emission into these modes with rate 2T. 



Hamiltonian must have the form 



We shall use the interaction picture to simplify the calculations in chapter El 

3.7 Decoherence and dissipation - the "quantum 
jump" approach 

In a real experiment, the cavity will not be perfect, photons in the cavity will escape 
by leakage through the mirrors, and the atom will spontaneously emit photons out 
of the cavity due to the coupling of the atom to the modes of the electro-magnetic 
field of the free space on either side, as illustrated in figure 13.21 

The dynamics of the atom-cavity system under these conditions can be described 
by a Master equation. This is an evolution equation for the reduced density operator 
for the atom-cavity system. The dynamics of the system under these forms of 
decoherence are Markovian, i.e., the environment does not retain a memory of the 
decay events. This means that the Master equation can be written in Lindblad 



-f^i.p. 



U\t)HU{t) +ih 



dWjt) 
dt 



U 



(3.24) 
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form [Lindblad76] , 



[HesP - pHlg] + JipJi 



t 



(3.25) 



P 



h 



where Jj are "quantum jump" operators, and Hes is called an "effective Hamilto- 
nian" for the system whose Hamiltonian without dissipation and decoherence would 
be H^: 



This decomposition of the Master equation has a pleasingly intuitive physical 
interpretation. The loss of energy from the system is described in terms of quantum 
jumps, which describe the loss of a single quantum of energy from the system. Each 
operator Jj embodies one type of jump. The non-Hermitian effective Hamiltonian 
describes the continuous Schrodinger-like evolution under the condition that no jump 
takes place. Its non-Hermitian character means that the norm of the state under this 
evolution will decrease, and this corresponds exactly to the probability that no jump 
has occurred. Thus the system evolution can be thought of as an average over many 
trajectories, each of which consists of continuous Schrodinger-like evolution and 
random quantum jumps^. This "quantum jump approach" [ZoUerST, Dalibard92] 
(see [Plenio98b] and [Zoller97] for introductory reviews) allows one to numerically 
simulate the Master equation using Monte-Carlo methods. 

For a two-level atom in a cavity, there are two major sources of energy loss, 

illustrated in figure 13. 2j leakage of the photons from the cavity at rate 2k, and 

spontaneous emission at rate 27. The jump operator for the photon leakage is 

v/2/ta, and for the spontaneous emission is y/2rf\A){B\. Although the construction 

of the Master equation in this way seems purely phenomenological, it can, in fact be 

derived from the physical decay processes by a fuller treatment (see, for example, 

[MeystreQQ], 3rd edition). 

^It must be stressed however, that no physical reality can be ascribed to these trajectories. 
There are infinitely many other equally valid possible "unravellings" of the Master equation. 




(3.26) 
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3.8 Generation of entanglement in cavity QED 

Cavity QED has proven itself to be an excellent testing ground for the validity of 
quantum mechanics with a high degree of precision. Many beautiful experiments 
have been carried out, and in recent years, entangled states have been created and 
verified. There are two principle approaches for the creation of entanglement in these 
systems. The first is a precise control of the coherent dynamics, such that the desired 
entangled state evolves through Rabi oscillation. This was the first approach to be 
implemented experimentally. Rauschenbeutel and co-workers generated Bell states 
between two modes in a single cavity using a Rydberg atom [RauschenbeutelOl] 
coherently interacting with each mode in turn. Their scheme could be adapted in a 
straightforward way to generate such states in spatially separated cavity modes. The 
same Paris group, which is led by Serge Haroche, have performed several further such 
experiments summarised in the following review articles [RaimondOl, HarocheOS] 
including the creation of an entangled state between two atoms and a cavity. 

In the second approach, measurement plays a key role. The first proposals 
of this kind were by Cabrillo and co-workers [Cabrillo99] and Plenio et al [Ple- 
nio99]. They showed that Bell states could be generated between atoms condition- 
ally, by driving them with a weak laser pulse, and subsequently detecting a photon 
spontaneously emitted by one of the atoms. The key to their scheme is that the 
measurement cannot in principle reveal from which cavity the photon originated, 
and the measurement thus projects the atoms onto the desired entangled state. 
Many refinements of this idea have been proposed, including the scheme by Bose et 
al. [Bose99] which is more relevant to cavity QED experiments will be discussed in 
detail in chapter Several enhanced versions of these schemes have been proposed 
recently, [Protsenko02,Feng03,Simon03,Duan03], which aim to combat some of their 
inherent drawbacks, namely that the schemes typically have a finite success proba- 
bility, often around 50%, and often do not work well when realistic imperfections, 
such as spontaneous emission and inefficient detection, are taken into account. It has 
recently been pointed out [Browne04] that such entangling operations can be used 
to generate entangled multi-qubit states known as cluster states [BriegelOl], which 
are a resource for universal quantum computation via single qubit measurements 
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alone [RaussendorfOS] . 

In the following two chapters, we describe two novel schemes for the generation 
of maximally entangled two qubit states between spatially separated cavities. Both 
of these schemes fall into the second of the above categories, and rely on the outcome 
of a measurement for the generation of the desired entangled state. In chapter |31 
we introduce a method of generating a Bell state between modes of two spatially 
separated cavities, whose robustness against experimental error stems from the fact 
that it relies on very short interactions between the cavities and a mediating atom. 

In chapter El this principle is taken a step further, where a single brief interaction 
is replaced by continuous weak driving, to produce a scheme for entangling two ions 
or atoms trapped in spatially separated cavities. This can, in contrast to previous 
schemes, achieve success probabilities arbitrarily close to 100%, and is robust against 
experimental imperfections. 
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Chapter 4 



Generation of Entanglement 
between two Optical Cavities 



4.1 Introduction 



In this chapter, we describe a novel method [BrowneOSb] of generating entangle- 
ment between two spatially separated optical cavities. The goal is the creation of a 
maximally entangled Bell state in the photon numbers of a mode in two spatially 
separated cavities. Let us label the two cavities A and B and label the Fock states 
of the chosen mode \n)j_ and respectively. The state which is to be generated 
can thus be written 



The scheme can be summarised as follows. A beam of atoms is employed, which 
have a transition resonant with modes in two identical cavities. The two cavities 
are positioned as in figure 14.11 such that their axes are parallel, and their centres in 
alignment. The cavities must be prepared in their vacuum state. Due to the high 
energies of their photons, cavities with optical resonant frequencies, known as optical 
cavities, are naturally very close to this state at room temperature. In the case of 
microwave cavities, on the other hand, this requires cooling to low temperatures. 

The first atom of the beam is prepared in the excited state of the resonant tran- 




(4.1) 



58 



4.2 Effective interaction time 



Atom Oven Collimator Detector 

State ^ Cavity A Cavity B 

Preparation 




Atomic Path 




Figure 4.1: A schematic diagram of the layout of the scheme. 

sition, which we denote as above, and then passes along a path through the 
centres of each cavity. Upon leaving the second cavity the atom's state is immedi- 
ately measured. If the atom is now in state 1^4), the lower state of the transition, 
and if the conditions explained below are satisfied, the entangled state l^"*") has been 
generated in the cavities. If the atom is measured to have remained in state \B), 
the procedure must be repeated, with a further atom. Let us examine the scheme 
in more detail, and explain how it works. A useful concept in the description of the 
atom-cavity systems after the atom's passage is the total effective interaction time 
between the atom and each cavity. 

4.2 Effective interaction time 

Since the interaction Hamiltonian of an atom interacting with a cavity 

H,nt = g{v)a\B){A\ + h.c. (4.2) 

is dependent on the position r of the atom, when the atom is in motion through 
the cavity this Hamiltonian becomes time- dependent. The Schrodinger equation 
for a time- dependent Hamiltonian can often not be solved analytically, and this is 
especially true in a situation like this, where the time-dependence is a functional 
of a yet unspecified path r(t). Nevertheless, in the special case that the detuning 
between the atomic transition and the cavity is zero, the Schrodinger equation can 
still be integrated analytically since the dressed states of the system do not depend on 
\g{j)\ and thus remain time-independent (see equations ()3.15|) and fj3.16|) ). Since 
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the mode-functions of the standing cavity-modes we are considering are real, the 
phase of g{v) remains constant along its path. Thus, for simplicity, and without loss 
of generality, we assume that g{v) is real. 

If we look, for example at the equation of motion for the system initially in the 
product state we see that is has a similar form to equation ()3.2()j) 

|^(t)) = e-^'^("+i)* (-1 sin (^f^^+l [ \g{t')\dA \A,n+l) 

\ \ Jo J (43^ 



except that the simple product Qnt = ^/n + Igt is replaced by an integral over 
time, ^/n + 1 Jl^\g{t')\dt'. That this equation of motion satisfies the appropriate 
Schrodinger equation can be verified easily by differentiation. If we define the ef- 
fective interaction time T{t) = g{t')dt' / go, where go = g{r = 0) is the coupling 
strength when the atom is in the centre of the cavity, then the equation of motion 

= e— (-^ sin (^or(t)) \A, n + 1) + cos (^or(t)) \B, n)] (4.4) 



is now equivalent to that of a stationary atom, sitting in the centre of the cavity 
for duration T{t). By calculating the equation of motion for a system initially in 
state \A,n + 1) as well, one sees that the unitary evolution operator for an atom 
interacting resonantly with a cavity with effective interaction time r is 

U = exp[-igoT{d\B){A\+d^A){B\)] . (4.5) 

If an atom passes through the cavity on a straight path from some point much 
further away from the centre of the cavity than waist of the cavity mode, to some far 
point on the other side of the cavity, one can characterise the complete interaction 
by the total effective interaction time. 
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4.3 Outline of the scheme under ideal conditions 

Using the effective interaction time, we can now show very simply how the scheme 
works. On its path through both cavities the atom interacts resonantly with each. 
Let us label the total effective interaction times with each cavity r4 and tq respec- 
tively. The evolution of quantum state of the system after the two atom-cavity 
interactions can be written 



where we assume that the coupling strength go is the same in each cavity and 
where we label mode operators for each cavity a_4 and respectively. 

In the scheme we propose, the cavities are aligned as shown in figure l^!T] and the 
atom passes along the straight line through the centre of both cavities at constant 
velocity. This means that the effective interaction times between the atom and each 
cavity will be equal. It may, of course, be difficult to control the path of the atom 
with sufficient accuracy that the interaction times are exactly equal, and the effect 
of this is discussed in section 14.4.11 For now, however, let us assume that r4 and rg 
are both equal to some value r. In the limit when the effective interaction times are 
very small, i.e. when qqt <^ 1, Uab can be expanded to the first order in g^r, and 
takes the following form 



Thus, the state of the system after the atom has left the second cavity, {tp) = 
UAsli'mit), can be written in this limit 



When the state of the atom is now measured, the cavity modes are projected 
into one of two states, depending on the measurement outcome. If \B) is detected, 
the cavity returns to its initial state. This is important for a non-deterministic 
process, because it means that it can be repeated immediately from the same starting 



U = exp [-igoTBiaBlB) {A\ + h.c.)] exp [~igoTys,{d_A\B){A\ + h.c.)] 



(4.6) 



U^l-igoT {aj, + aB)\B){A\ + {a\ + al)\A){B\ 



(4.7) 




(4.8) 
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conditions. If the ground state \A) is detected, the cavity modes are now in the (un- 
normahsed) state {a]^ + ae)|'0cav)- Thus, if the cavities are initially in the vacuum 
state, the state generated in the cavities would be {a\ + a^)|0)^|0)i3 = |1)^|0)b + 
|0)yfi|l)B, which when normalised is the Bell state |^+) introduced above. 

If, following a successful run, one were to repeat the scheme immediately and 
carry on until n atoms had been detected in the ground state, in the limit that QqT is 
small, the state generated would have the form (a^ + ag)"|'0cav)- More specifically, if 
the cavities are initially in their vacuum state, the state would be (a^+ag)"'|0)^|0)B. 
This is equivalent to the state produced when a n-photon Fock state and a vacuum 
state are incident together on a 50:50 beam splitter. However, the numerical simu- 
lations we have undertaken have suggested that the fidelity of states produced via 
this method would decrease swiftly with increasing n. This is due to two reasons. 
Firstly, the short interaction time approximation becomes worse when higher photon 
numbers are present in the cavities since the time scale of the interactions is faster 
(the Rabi frequency scales with ^/n+ 1). Secondly, when more than one photon is 
in the cavities, a failed attempt does not reset the state of the cavities to the state 
before the run, so, as more repetitions are made, the fidelity of the final state gets 
worse and worse. For these reasons, this does not appear to be a good scheme for 
the generation of such states. 

However, as we will show below, the generation of single-photon Bell states is not 
affected by these problems. Firstly, we find that the fidelity of the generated state 
remains close to unity for values of QqT much greater than the above approximation 
would be valid. Secondly, "failures" reset the state of the cavities to the vacuum 
state exactly, so the fidehty of the state generated is unaffected by the number of 
runs which were required. In the limit that QqT tends to zero, the probability of the 
detector measuring a ground state is approximately 2{goT)^, which, in this limit, is 
vanishingly small. In a practical scheme, one would need to work in a parameter 
range where the probability of success was high enough that few repetitions would 
be required to achieve a successful run. However, this is the case for higher values 
of goT, where the above short-time approximation would no longer be valid. For- 
tunately, starting with the simple pure state |i?)|0)^|0)B, it is straightforward to 
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solve the Schrodinger equation and calculate exactly the state of the system after 
the interactions have taken place. This state lip), under the assumption that the 
interaction times r are exactly equal is, 

1^) =cos'igoT)\B)\0)A\0)t3-icos{goT)sinigoT)\A)\0)A\l)B , , 

(4.9) 

-isin((7or)|A)|l)^|0)B. 

Let us consider a measurement of the atom's state. If the excited state is detected 
the state of the cavities is projected back to the vacuum state, independent of r, 
as in the approximate case. If the ground state is detected, the following entangled 
state is generated in the cavity, 

l^cav) = cos((?or)|0)^|l)B + |l)^|0)e, (4.10) 

where normalisation has been omitted. As g^r approaches zero, this tends to the 
desired state I^E'"*'). One can quantify how close this state generated is to I^E'"^) in 
terms of the fidelity F = |(^+|^cav)P, 

2 cos''((7ot) + 1 16 

Note that there is no term proportional to {gor)^ since the series expansions of 
cos{goT) and {cos'^{goT) + 1) are proportional up to this power of (gor). The expres- 
sion is plotted in figure We see that the fidelity remains very close to unity for a 
surprisingly large range of g^r. For example, the fidelity remains above 1 — 4 x 10~^ 
for g^T = 0.5. Additionally, the function is extremely flat in the wide range of val- 
ues between goT = and goT = 0.5. Therefore the scheme would be insensitive to 
variations in the interaction times within this range. 

For the scheme to be useful, the probability of a single run leading to the suc- 
cessful generation of an entangled state, Psuccess, needs to be high enough that pro- 
hibitively many repetitions are not required. 

^success = 1 - cos\goT) = 2{g^Tf + O {{gorf) (4.12) 
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Figure 4.2: The probability that a single run will lead to the detection of the atom in 
its ground state, and thus the successful generation of a Bell state, plotted for values of 
goT from to 1. 
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Figure 4.3: The fidelity of the Bell state generated in the cavity modes after a successful 
run, plotted against qqt with values from to 1. Note that the values of goT corresponding 
to a full Rabi oscillation in each cavity is goT = vr. 
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The success probability is plotted in figure 14.21 For qqt = 0.5, this probability 
is approximately 0.4, so a successful run would probably be achieved in 2 or 3 
repetitions. The optimal parameter range for the scheme depends upon the fidelity 
of the state required. The higher the value of qqt chosen the higher the success 
probability will be but the lower the fidelity of the state generated. If, for example, 
fidelities of 0.95 were acceptable, and a minimum success probability of 0.5 were 
desired, the scheme could operate between approximately g^r = 0.6 and goT = 0.9. 
The exact value of goT, however, can lie anywhere in this range, so no fine tuning of 
the interaction time is required. Note that even if extremely high fidelities such as 
0.999 are required, the necessary parameters (up to g^r = 0.35), still allow a success 
probabihty of up to 0.22, meaning that a successful run would probably be reached 
after 4 or 5 repetitions. 

4.4 Practical considerations 

The above description is quite general and no particular type of atom or cavity 
has been specified. A number of aspects must be considered in the choice of a 
physical system to implement the scheme. Firstly, it would be desirable that the 
entangled states, once generated, would be as long-lived as possible. This favours 
microwave cavities over optical cavities, since the lifetime of a photon in an optical 
cavity is currently at the very most one microsecond [Lange02], whereas microwave 
cavities with a photon lifetime of a millisecond have been made [RaimondOl]. It 
takes an atom travelling at 500ms~^ around 20/is to traverse a typical microwave 
cavity, which is much smaller than this lifetime, so there would be sufficient time for 
further atoms to probe and interact with the cavity mode before the entangled state 
has dissipated. If a microwave cavity is used, an atom with a microwave transition 
is then needed. Microwave transitions occur in atomic fine and hyperfine structure, 
however dipole transitions between these states are forbidden, and their interaction 
with the cavity mode would be much too weak to implement this scheme - the 
interaction times that would be required would be much greater than the photon 
lifetime of the cavity. Rydberg atoms, on the other hand, although more difficult 
to prepare, have large dipole moments and thus would interact strongly with cavity 



65 



4.4 Practical considerations 



modes. 

Experiments have been carried out, for example by the Paris group of Haroche 
and coworkers [RaimondOl], which have parameters close to that required in this 
scheme. Recall that our scheme requires that the product of parameters go and r is 
at the minimum 0.2 and maximally 0.5 to 0.8, depending on the fidelity of Bell state 
one wants to generate. In the Paris experiments, a Rydberg atom interacts with a 
microwave cavity Rydberg atoms interact resonantly with microwave cavity modes 
with a Rabi frequency of approximately 47kHz. This means that go = 470007rs~^ = 
1.48 X lO^s^^. Atoms in an atomic beam from an oven source travel at speeds of 
the order of hundreds of meters per second. In the Paris experiment atoms with a 
speed of 500ms~^ are selected. This means that the effective interaction time is such 
that a single Rabi oscillation is performed, i.e. goT = tt. This is a factor of 4 to 8 
lower than the parameter range for our scheme. It would be difficult to lower r by 
using faster atoms, since the velocity of the atoms scales with the square root of the 
atom oven temperature, so a lower go would be required. This could be obtained be 
using a larger cavity, and since go scales with 1/VV = so a cavity of mirror 

separation 3 or 4 times as great as in the Paris experiments would lead to a effective 
interaction times in the required range. Thus Rydberg atoms and microwave cavities 
could be employed to implement the scheme. 

The disadvantage in using Rydberg atoms is that, at present, the efficiency of 
state detection schemes is much lower than unity. In [RaimondOl], for example, they 
report a detection efficiency of 40%. We will discuss the implications this has for 
the scheme in subsection 14.4.21 First, however, we consider the effect on the fidelity 
of the entangled states produced in the scheme if the path of the atom through the 
cavities is not well controlled and deviates from the line through the centre of the 
cavities. 

4.4.1 The atomic path 

If the effective interaction times with both cavities are not exactly the same, this 
can reduce the fidelity of the entangled state produced by a successful run of the 
scheme. In our discussion above, we assumed that the two interaction times were 
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exactly equal. In practice, however, it could be difficult to control the path of the 
atom so precisely. Let us consider ffist the effect that differing interaction times 
would have in the fidelity of the Bell state generated by a successful run of the 
scheme. Let the effective interaction between the atom and cavity A and the atom 
and cavity B he t and r(l — e), respectively. We can rewrite equation ()4.11|1 . to 
take the differing interaction times into account, and find the following expression 
for the fidelity, 

^ ^ 1 ^ cos(^or) sin(^or) sin(^or(l - e)) 
2 cos2(5(or) sin^ (fi'or(l - e)) + sin^igor) 

The fidelity is plotted for qqt = 0.5 as a function of e in figure ^31 The asymme- 
try of the plot is partly an artifact of the choice of parameterisation, but, if we take 
this into account, by plotting F against ln(l — e), as in figure l^75| we see that the 
asymmetry remains. This is due to the asymmetry in the scheme itself regarding the 
interactions with the two cavities. When the atom enters the ffist cavity it is always 
in the product state |-B), whereas, when entering the second it is always entangled 
with the first cavity. This makes the interaction with cavity B slightly weaker, and 
is the reason that the maximal value of F occurs when e has a small negative value. 
The slightly longer interaction time compensates for the interaction being slightly 
weaker. 

Since it is reasonably easy to position the cavities in the desired place to very 
high precision, let us assume that they are perfectly aligned in the layout illustrated 
in figure W7I\ The factor which would be harder to control would be the atom's path 
through the apparatus, since it would be travelling on a ballistic path after being 
ejected from a heat oven, and although such atomic beams can be highly coUimated, 
there will generally be a small amount of spread in the transverse direction, meaning 
that the atom's path will diverge slightly from the central path. 

To calculate how much this would affect the interaction times we consider the 
cavity geometry. We assume that the cavities are in a Fabry-Perot configuration 
as discussed in section 13.11 Since the atoms will be assumed to remain close to the 
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Figure 4.4: The fidelity of the Bell state generated in the cavity, if the effective inter- 
action times are r and r(l — e) for the interactions with cavities A and B respectively, is 
plotted here ctS cl function of e for 50''" = 0-5. 
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Figure 4.5: The fidelity of the Bell state generated in the cavity plotted again as a 
function of e for gQT = 0.5, this time plotted against ln(l — e). 



2;-axis of each cavity, we may use the approximate mode function 



u{x,y,z) 



2tiz 



cos 



(4.14) 
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Figure 4.6: A general straight line path through the cavities can be defined in terms 
of the parameters yo, zq, (p and 9, where distances are measured in meters and angles 
in radians. This figure illustrates zq and 0, the other two parameters, yo and (/>, are 
equivalently defined in the perpendicular x — y-plane. The distance between the exit of 
the collimator and the centre of cavity A is Dq. The distance between the centres of the 
cavities is Di. 

where wq is the mode waist. The z-axis hes along the hne connecting the centres of 
the two mirrors and the origin is in the centre of the cavity. The waist of a Gaussian 
mode, Wo, is a function of the cavity geometry and the field wavelength. 



2 \JL{2R-L) 
< = ^^A^ (4-15) 

where, as above, L is the separation between and R the radius of curvature of the 
mirrors. If the atom travels along the central axis of the cavity with constant speed 
f , the effective interaction time is ^/^^wo/v. Let us consider a general straight path 
through the system, which can be defined in terms of four parameters, yo and zq, 
the initial displacements from the central line in the y and z directions, and (f) and 
9 the angles between the atomic path and the central line in the y and z directions, 
illustrated for zq and 9 in figure W7\ We can calculate the effective interaction times 
between the atom and cavity as it travels along this path with speed v and find, 
that the interaction time with each cavity, in terms of displacements 5y and 5^ which 
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differ for each cavity, is 



Teff = 



f cos6' 



COS k6z — k6,, 



exp — 




4 



tan' 



2^1 



(4.16) 



For tlie interaction with cavity A, 6y = ?/o + cos(0)Do and 6z = zq + cos{9)Dq, and for 
the interaction with cavity B 6y = yo+cos{(f)){Do+Di) and 6z = zo+cos{9){Do+Di). 
We can use these expressions to calculate e in terms of these parameters, written 
here to the second order in i/q, zq, and 9 and their products. 



Let us discuss the constraints this would have on the coUimation of a typical ex- 
periment. In the cavity QED experiments of Haroche and coworkers [Brune96, Rai- 
mondOl], Rydberg atoms interact resonantly with microwave cavities. In a typical 
experiment, a cavity mode with waist wq = 5.97mm and wavelength A = 5.87mm 
is employed. For e to be small, the quantities in parentheses in equation ()4.17p 
must be much smaller than wq and \/^/27i = 1.32mm. This means that the atomic 
beam must be coUimated so that the effective beam radius is much smaller than 
this distance. In [Brune96] an effective beam radius of 0.25mm is reported. If we 
assume from this that, in the worst case, this would mean that yo, Zq ^0.25mm and 
Di(j), Di9 ^0.25mm, we can estimate that e would be less than 0.2. If qqt = 0.8 and 
e = 0.2, this would corresponds to a reduction of the fidelity of the entangled state 
produced by a successful run from 0.96 if both interaction times are exactly equal 
to 0.93. Equation ()4.17|) also provides another reason why optical cavities would 
be unsuitable for the scheme. The typical waist of an optical cavity tends to be 
much smaller than that of a microwave cavity, so the demands on the atomic beam 
coUimation required if optical cavities were used would be extremely high. 



e 



\ [{D,<f)f + (Z}i0) (2/^00) + (2yo)(A0)] 



+ ^ [{D,9f + {D,9){2Do9) + {2zoKD,9)] . 



(4.17) 
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4.4.2 Detector efficiency 

Our analysis in the previous section assumes that the atomic state detector has per- 
fect efficiency. In practice, the detection efficiency will be less than unity. Indeed, 
as mentioned above, current state detection methods for Rydberg atoms have an 
efficiency of just 40% [RaimondOl]. The detection process ionises the atom, de- 
stroying the state, so increased efficiency cannot be obtained by placing detectors 
in series. A single detection failure will disrupt the scheme, since it will cause a 
mixed state to be created in the cavities. The scheme must then be halted, and one 
would then have to wait until this state would have dissipated from the cavities, 
and the cavities have returned to the vacuum state. Otherwise, if further atoms are 
sent through the cavities immediately, they will interact with the mixed state, and 
any "successful" run, will have generated a mixed state with much reduced fidelity. 
Rather than halting the flow of atoms through the cavities, their interaction with 
the cavities could be prevented for the cavity dissipation time by the application of 
an electric held to the system, to create a Stark shift in the atoms such that they are 
no longer resonant with the cavities. Therefore, one would like a detection efficiency 
high enough that the probability of a detection failure, during the typical number 
of runs needed before a successful ground state measurement is made, is low. The 
mean number of runs to generate the entangled state in the cavity, if the detector 
were ideal, is simply the inverse of the success probability. The probability that the 
detector works every time during the process is therefore simply Di/^su«ess^ which 
for all values of Psuccess is less than or equal to D. Therefore, for a rehable scheme, 
a high detection efficiency would be desirable. 

For example, let us consider an implementation of this scheme with Rydberg 
atoms travelling through the cavities such that QqT = 0.5. The success probability 
-Psuccess is 40.7% and with current detectors with efficiency 40%, Pdct would be around 
10%. This would mean, typically one would have to repeat the whole process, 
including preparation of the cavities, 10 times before it could reach its successful 
conclusion. However, in light of the fact that the cavity dissipation time is of the 
order of milliseconds, and the time taken for each run much less than this, even in 
this case, the time needed to repeat the scheme enough times to generate the Bell 
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state would be a fraction of a second. 

4.5 Conclusion 

We have described a scheme for the generation of high-fidehty Bell states between 
two spatially separated cavity modes. The scheme is non-deterministic, but we have 
shown that within the range of parameters goT between approximately 0.3 and 0.9, 
fidelities higher than 0.95 are obtained, with success probabilities for a single run 
greater than 1/5 for this entire range. 

The most appropriate physical system to implement this scheme would be a 
combination of microwave cavities and Rydberg atoms. The low detection efficiency 
for Rydberg states would increase the number of times which the scheme would need 
to be repeated before a successful run, and require extra time after each detection 
failure to allow the mixed state produced in the cavities to dissipate. Nevertheless, 
the scheme would still be successful within a reasonable number of repetitions. The 
scheme requires that the atomic beam used is highly coUimated otherwise the fidelity 
of the states produced may be degraded, but the coUimation which has already been 
achieved in current experiments is high enough that this effect would be small. 

There are some disadvantages in generating entangled states in cavity modes. 
Firstly, they have a short life time which means that either the state would have to 
be utilised within a short time-window, or it must be mapped onto a fresh atom. 
This mapping can be achieved, however in the same experimental framework, except 
that the atom must be introduced in the lower state 1^4) and remain in the cavity 
for an effective Rabi oscillation. An alternative approach would be to work with 
atoms or ions trapped within the cavities via external fields. In the next chapter we 
consider such a situation, and demonstrate again the advantages that an approach 
based on weak interactions combined with measurements can have. 
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Chapter 5 

Robust Creation of Entanglement 
between Ions in Spatially Separate 
Cavities 

5.1 Introduction 

In the previous chapter, we introduced a scheme for entanglement generation, where 
high fidehties are obtained, because we work in a regime where the interactions in- 
volved are very brief. However, an unavoidable side-effect of the use of this technique 
is that the probability of success for a single run is low. In this chapter, we use similar 
ideas to greatly enhance a previous proposal for the generation of a maximally entan- 
gled two-qubit state between ions trapped two spatially separated cavities. Various 
schemes for trapping ions and atoms have been proposed and experimentally im- 
plemented [Paul90, Liebfried03] including trapping inside cavities [Ye99, KellerOS]. 
Our proposal uses once again the principle of using a slight interaction plus a mea- 
surement to achieve experimental robustness. This time, however, by applying a 
driving field continuously, but weakly, this robust scheme also exhibits arbitrarily 
high success probabilities. 

The structure of this chapter is as follows. First we will briefly summarise the 
proposal of Sougato Bose and co-authors [BoseQQ] and discuss its advantages and 
drawbacks. Then we will introduce and discuss our improved version of the scheme 
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Figure 5.1: Internal level scheme of the ions. Levels \A) and \C) are coupled via a non- 
resonantly interaction with a cavity mode. Levels \B) and \C) are coupled by an external 
"classical" driving field. Spontaneous decay is indicated by arrows from \C) to \A) and 
\B). 

[BrowneOSc] and demonstrate it's robustness against experimental imperfections. 

5.2 System Hamiltonian 

Before describing this proposal in detail we first derive a Hamiltonian which describes 
the evolution of the the ion-cavity system, which is also applicable to our modified 
scheme. 

We consider atoms (or ions) with three energy levels \A), \B) and \C), sitting in 
a cavity such that a single cavity mode of frequency Uc couples with the transition 
between levels 1^4) and \C) with detuning A. The transition between levels \B) and 
\C) is driven by an external field with frequency ui which we can treat classically. 
Neglecting decay processes, we can write down a full Hamiltonian from the system, 
using the results from chapter El 

H =hA\C){C\ - nuJc\A){A\ - hui\B){B\ + Uu^a^a 

n (5-1) 

+ h{g\C){A\a + h.c.) + h-{e-''^''\C){B\ + h.c.) , 
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where g is the atom-cavity mode couphng constant, and is the classical Rabi 
frequency of the interaction with the external laser field. Let us simplify this by 
transforming to an interaction picture with f^int = hA\C){C\ +h(^g\C){A\a + h.c.) + 
h{Q/2){e-''^^'\C){B\ + h.c.) 

Q (5-2) 

= hA\C){C\ + n{g\C){A\a + h.c.) + h- {\C){B\ + h.c.) . 

The schemes below operate in the high detuning regime, where A is much greater 
than the Rabi frequencies of each interaction. We may therefore simplify the Hamil- 
tonian further by adiabatic elimination of level \C). To perform this, first we diag- 
onalise the Hamiltonian, as we did in chapter El for the two-level system. We write 
Hamiltonian is rewritten as a sum of operators Hn each acting on the sub-space 
spanned by \A,n + 1), \B,n) and \C,n), and write Hn in matrix form as 



Hn = h 



/ g*VnTT \ 

n/2 

\ gVn + 1 ^1/2 A J 



(5.3) 



In the limit that A ^ g,^, the eigenvalues of this matrix become approximately 
the following values, Ei = 0, E2 = h{4:\g\'^{n + I) + il^)/(4A) and E3 ^ A, with 
corresponding eigenvectors 



= \/ 4|,|> + i) + ff (-"l-'-" + ^> + ^^^l^-">) • f^'^' 



" + 1) + {^s'V^MAn + 1) + n\B, n)) , (5.5) 

\i'3}~\C,n). (5.6) 

Thus, in this limit, \C,n) is an approximate eigenstate of the full Hamiltonian, 
and does not couple to the n+ 1), \B,n)} subspace. Since the initial state 
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of the system in our scheme will be the in latter sub-space, we therefore simplify 
the Hamiltonian by restricting it to this sub-space, adiabatically eliminating \C,n). 
Casting the restricted diagonalised Hamiltonian back in terms of the product state 
basis vectors, we obtain 



r^Mgl (n + i) + n\ , , . , 



4A 



n 



n\3l\A){A\a^a + h^\B){B\ + fh^\B){A\a + h.c. 



(5.7) 



Note that if system is prepared in the {\A, 0)} sub-space and the classical 

field strength is chosen such that Q = 2g then this Hamiltonian is exactly like a 
resonant two-level Jaynes-Cummings Hamiltonian with Rabi Frequency 2g'^/A. 



5.3 An earlier proposal 

Here, we introduce a proposal [Bose99] for the generation of a maximally entangled 
state between two ions trapped in spatially separated cavities, on which our improved 
scheme is based. Although this proposal was primarily presented as a teleportation 
scheme by Bose and co-workers, we focus here on the author's modified version 
in which a maximally entangled state is prepared between two ions (or atoms) 
trapped in two spatially separated cavities. Bose and co-authors considered the 
setup illustrated in figure 15.21 Two cavities each contain a trapped ion (or atom) 
driven by an external field, as described in figure 15.11 and in the previous section. 
One mirror in each cavity is semi-silvered such that photons leak out at rate 2k. 

Each atom is prepared in state \B) and the cavities are cooled to their vacuum 
state, before the classical driving fields are applied. Driving fields with strength 
Q = 2g are applied for exactly one quarter a Rabi cycle, that is for duration n/Ag'^, 
and then the driving field is switched off. This drives each ion-cavity sub-system 
through a quarter of a Rabi oscillation into the state, 

\i^) = \ll{\Bm-z\A)\l)) (5.8) 
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Figure 5.2: We consider a set-up in which individual ions are trapped inside two spatiahy 
separated optical cavities A and B where they are also driven by an external classical 
field. Photons can leak out of the cavities and are then mixed on a beam splitter B.S. and 
subsequently detected by photo-detectors Di and D2. 

where |0) and |1) are the zero and one photon states of the cavity field. Adding 
subscripts A and B to denote the cavity, the combined state of both ion-cavity 
systems is 

\^) = U\B)^\B)s\0)a\0)s-\A)a\A)s\1)a\1)b 

^ (5.9) 

-t[\B)A\A)s\0)A\l)B+\A)A\B)B\l)A\0)B]] ■ 

For simplicity here we assume that the leakage rate is much slower than the Rabi 
frequency of the state preparation step above, and neglect it (a fuller treatment may 
be found in [Bose99]). 

One can model this mathematically using a Master equation in the quantum 
jump approach, as described in section 13. 7[ however, we can explain Bose et al's 
scheme qualitatively without resorting to this level of analysis. Let us simply assume 
that the photons will leak out of the cavity at some time, where they will mix on 
a beam splitter and be detected by detectors Dl and D2. Since the cavities leak 
at the same rate, when a single photon is detected, one cannot distinguish from 
which cavity it originated, and this lies at the heart of the scheme's ability to create 
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entanglement. 

The action of the 50:50 beam sphtter on the four possible input basis states 
|0)|0), |0)|1), |1)|0) and |1)|1) is the following transformation 

|0)|0) 
|0)|1) 
|1)|0) 

One waits for a length of time such that it is very likely that all photons in 
the cavities have leaked out and, if present, been detected. One can interpret the 
detection outcome as a measurement on the initial state of the cavity. If one and 
only one photon is detected, then this is equivalent to having performed a joint 
projection on the states of the two cavities onto |0)^|1)b + ^|1)^|0)H or i|0)^|l)B + 
|l)yi|0)B depending on which detector clicked. This projection leaves the two ions 
in a maximally entangled state — i\B)yi^\A)is or —i\A)j[\B)j^ + \B)^\A)q. 

Neglecting imperfections in the state preparation due to the cavity leakage, and 
assuming perfectly efficient detection, the success probability is 50%. Half the time 
either zero or two photons will be detected, which leave the ions in a product state. 
The scheme must then be repeated. Note that a detector which can discriminate 
between different one or two photons is needed here, a great technological challenge. 

In [Bose99], the authors simulate the performance of the scheme with a full 
treatment of the cavity decay and consider inefficient detectors. They find, in their 
more thorough treatment, that the success probability is lowered by the small chance 
of photon leakage during the preparation stage. They also find that low detection 
efficiencies are particularly damaging, as they mean that states which should have 
led to two photons being detected generating product states are wrongly identified 
as successes, which badly degrades the fidehty and the degree of entanglement of 
the states generated. 

Using the same basic setup as this scheme, we now use the principle of employing 



H0)|0) 



(|0)|1) + ^|1>|0)) 
(^|0)|1) + |1)|0)) 
(|0)|2) + |2)|0)). 



(5.10) 
(5.11) 

(5.12) 

(5.13) 
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weak or brief interactions and measurement to give a new proposal [BrowneOSc] 
which, (i) succeeds with 100% probabihty under ideal conditions, (ii) allows the 
achievement of high fidelity entanglement outside the strong coupling regime upon 
the detection of a photon, (iii) is robust against detector inefficiencies and absorption 
losses in the cavity mirrors and (iv) can be adapted, with the same efficiency, to 
entangle ions trapped in a single optical cavity. 

5.4 Weakly driven approach 

We begin the discussion of our new scheme by first analysing in more detail the 
evolution of each of the driven ion-cavity sub-systems including cavity leakage. Our 
initial analytic treatment will be valid for the so-called "strong-coupling regime" 
where spontaneous emission from the ions is neglected. After this we will consider a 
more experimentally realistic "weak-coupling regime" with significant spontaneous 
emission and also include finite detection efficiency in our model. 

First let us concentrate on a single ion-cavity sub-system. We assume here that 
light leaks through one of the cavity mirrors, with rate 2k. Adopting a quantum- 
jump approach, we write the Master equation for the ion-cavity system, 

p = ~i{HcsP - pHls) + 2Kapa^ (5.14) 
where we have defined the non-Hermitian effective Hamiltonian 

^eff = Had - ii^a^a (5.15) 

where Had is defined above in equation ()5.7|) . 

In contrast to [Bose99] , our proposal operates in what we call the "weak driving 
limit" such that the condition |^ ^ k is satisfied. Intuitively this implies that the 
rate of transitions between levels 1^4) and \B) of the ions will be slower than the rate 
of cavity decay. Since our scheme is initiated with both ion-cavity sub-systems in 
the product state |-B,0), this implies that the population in level \A) of the atoms 
will be small, unless a photon is detected. More precisely, in the weak driving limit, 
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under the condition that no photon is detected, the state of the sub-system will 
quickly approach the following state 

where x = — 

To see this, we diagonalise the non-Hermitian Hamihonian in equation ()5.15j) . If 
the sub-system begins with the ion in state \B) and both cavity and external modes 
in the vacuum state, it remains in the sub-space {|-B, 0), \A, 1)}. In the weak-driving 
limit we obtain the following eigenvalue, eigenvector pairs. 



El - ^<\x\') |ei) = y' 1) + \B, 0)) (5.17) 

E, ^M! _ le^) = ^-±-^(\A, 1) - x\B, 0)) (5.18) 

where x is defined above. 

Now, since this Hamiltonian is non-Hermitian, the eigenvalues are not all real 
and the eigenvectors are not orthogonal. Nevertheless, since the eigenvectors are still 
linearly independent, one can still represent an arbitrary state {ip) = a\A, 1)+P\B, 0) 
as a unique superposition of these (normalised) eigenstates^. Once decomposed in 
this form, the time evolution, under the condition that no photon leaks from the 
cavity, may be simply calculated as above by associating the factor e~*^'*/^ with 
each eigenstate, where Ei is the associated eigenvalue. Since Ei may be complex, as 
well as the usual oscillatory part of this factor, the imaginary part of Ei contributes 
to an exponential decay factor. 

Now, in the weak-driving limit of our non-Hermitian Hamiltonian, the imaginary 

part of El is tiny, but E2 has the large imaginary part —in. Thus the state will 

converge towards eigenvalue |ei) at rate 1/k, with a success probability dependent on 

the overlap between the starting state and |ei). Since, in the weak-coupling regime, 

^Note, however, the weights associated with each eigenvector are no longer simple overlaps with 
lip) but have a more complicated form. 
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the overlap between |ei) is close to unity, the norm of the conditional state remains 
very close to unity during this time scale, so the chance of this state preparation 
stage of the protocol faihng via the leakage of a photon is negligible. 

Thus, after a short period, both ion-cavity sub-systems have attained the state 
|ei). Wc now recall that the light escaping from the two cavities is mixed at a 
50:50 beam-splitter, and then passes into the detectors Dl and D2. Now that we 
are considering both atom-cavity subsystems, let us give operators the suffices A 
or B according to which subsystem they belong. Due to the linear transformation 
imposed by the beam splitter, it is impossible for the detectors to distinguish from 
which cavity the detected photon originated and a detection at detector Dl will 



now correspond to the jump operator Ji = ^/T/2{id_A + clb) and a detection at D2 
to J2 = v^l/2(a^ + ias). We write the Master equation for the combined state of 
both atom-cavity systems as follows. 



Ha,b P - pH\q + 2KaApa\ -\- 2KdBpal (5.19) 



where H^^s = {Hes)A + {Hes)B)- One can rewrite the jump part of the equation 
in a form more suitable to the detection scenario by using the identity 2KaApa\ -\- 
= 2kJipJ\ + 2kJ2pJ2- 
Let us assume that each sub-system begins in state |ei). The state of the system, 
after a very short time At will be 

p(Ai) = p(0) + At [-i \^Ha,b p(0) - p(0)^Xs] + 2«Jip(0)4 + 2«J2p(0)j|) 

(5.20) 

which in the weak coupling limit, where x <^ 1, gives to a close approximation, 

p(Ai) (1 - 4K|x|2Ai)|V'o)(V'o| + 2k\x\^M\'4)i){'4)i\ + 2K|a;|^At|V'2)(V'2| (5.21) 
where |^o) = |ei)^|ei)B and 

^ |5,0)^|A0)^ + ^|A,0)^|5,0)b 



x'^"' V2(l + |x|2) 



(5.22) 
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and likewise 



\'ip2) = — IV^o) 



+ x{\A,0)A\A,l)B + t\A,l)A\A,0)B 



t\B,0)^\A,0)s+\A,0)^\B,0)B 

(5.23) 



which are normahsed in the hmit that x is small, when both of these become very 
close to a maximally entangled state in the ions' internal degrees of freedom. 

Thus with probability 4/t|xpAt one of the detectors clicks, and a state very close 
to a maximally entangled state is created between the ions. Otherwise, there is no 
click and the state is reset to \ipo)- If one models the detectors as taking repeated 
measurements at some short time At apart, then in each time window the probability 
of a click is 4K|xpA)f: multiplied by the probability that no detection has occurred 
so far. Taking At to zero, the probability of the first click happening at time t is 
P{t) = 4k|xP exp[— 4K|xpt] and thus the mean time before the first detection event 
will be 

1 kA 



= X = = (5.24) 

In the weak driving limit, x is much smaller than 1 and can be made arbitrarily 
small by increasing the detuning and weakening the driving fields. Thus a maxi- 
mally entangled state of arbitrarily high fidelity can be generated with unit success 
probability. 

However, this result is still only valid in the strong coupling limit as we have 
so far neglected the effect spontaneous emission. We have also assumed perfectly 
efficient detectors which are not experimentally available. We can investigate how 
our scheme performs taking these effects into account by numerical simulation. 



5.4.1 Spontaneous emission 

Let us first consider the effect of spontaneous emission. Since spontaneous decay will 
principally occur from state \C), we can no longer use the adiabatically eliminated 
Hamiltonian, and return to the Hamiltonian in equation ()5.1|) for each cavity. We 
add spontaneous emission jump operators to the Master equation as described in 
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section EIZl assuming that two decay channels exist, from \C) to \A) with rate 7^ 
and from \C) to \B) with rate 7^. We can now use the quantum jump approach 
to calculate the evolution numerically using a Monte Carlo simulation. Since now 
sometimes the photons will be lost by spontaneous emission, sometimes neither 
detector will click. Thus, one must choose a "waiting time", T^, the length of 
time the detectors will be monitored before the attempt is judged to have failed, 
and the system re-initiated. This in turn, gives us a meaningful concept of success 
probability for the scheme. There are also two parameters relevant to our simulation 
technique. The first is the simulation time step. This is the length of time between 
the "coin tosses" which choose whether a quantum-jump has taken place. For a 
good simulation this time must be small compared to the decay times of the norms 
under the effective Hamiltonian evolution. We find that these remain on the order of 
^av = 1/4k|xP, even for large spontaneous decay rates, as expected and thus choose 
time steps of Tav/10^. Secondly, we need to repeat the simulations enough time 
such that a fair average is obtained. Clearly, there is a trade-off between accuracy of 
simulation and calculation time. Since we are more interested here in the qualitative 
behaviour of changing parameters than an exact quantitative calculation, we choose 
10^ repetitions, which gives sufficient accuracy and reasonable calculation times. 

The results of a simulation of the effect of spontaneous emission are plotted in 
figure Eini for the the parameters detailed in the figure caption. As expected, we see 
that finite spontaneous emission rates cause a linear reduction in the fidelity and 
the success probability. This corresponds to events where spontaneous decay has 
driven one or both of the cavities into the state |A, 0). This state is an eigenstate 
of the single cavity Hamiltonian and the system remains there. If this happens in a 
single cavity, any following detector click will lead to the undesired, and unheralded 
production of a product state, leading to a reduction in the fidelity of the state 
generated, and if this happens in both cavities no detection will occur, leading to 
the reduction in success probability. 

Working with a weak coupled system with strong spontaneous emission, one 
would expect that the fidelity of the generated state be maximised by reducing the 
waiting time since this reduces the chance that a spontaneous decay has occured. 
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Figure 5.3: The success probability (solid green line) and the deviation from unit fidelity, 
are plotted versus the spontaneous emission rate, which is assumed to be equal for both 
decay channels. The parameters chosen for the simulation were = 2g, k = Wg, A = 205. 
A time-step of O.lg and a waiting time of = 100 /g were used. The plot has been 
obtained from a quantum jump simulation where each point is the result of an average 
over 10^ runs. 

as a trade off against success probability. This confirmed by the simulation results 
plotted in figure EiH 

5.4.2 Detection efficiency 

Another important experimental imperfection is non-unit detection efficiency. In 
fact, not only are current photo-detectors unable to detect all incident photons, but 
other effects, such as absorption of photons in the cavity mirrors and the photons 
escaping from the cavity in the wrong direction will all contribute to give an effective 
detector efficiency rj much lower than one. There is a simple way to model such a 
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Figure 5.4: The success probability (solid green line) and the deviation from unit fidelity, 
are plotted versus the waiting time T^/TaV for VL = 2g, k, = 10(7,731 = 732 = O.li^ and 
A = 200(7. Each point is the result of an average over 10^ runs of the scheme. 

situation. The case of imperfect detectors is modelled as follows, the usual projection 
operator \n){n\ is replaced by 



Each term in the sum corresponds to {k — n) photons being missed or absorbed 
before the detector. 

It is straightforward to incorporate imperfect detectors into our Monte Carlo 
simulation. We expect that our scheme will exhibit robustness against such an 
imperfections since, in the weak-driving regime, there is only a very small probability 
that photons will leak from both cavities, which is the only event which will decrease 




k=n 



(5.25) 
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Figure 5.5: The success probability (solid line) and the deviation from unit fidelity, are 
plotted versus the detector efficiency for Q = 2g,K = IDS', 731 = 732 = O-lg and A = 20g 
and a waiting time of T = 100/ g. Each point is the result of an average over 10^ runs of 
the scheme. 

our state fidelity, by contributing an admixture of the product state |y4)|y4). Indeed, 
the simulation backs up this view. In Fig. 15.51 we have plotted both the success 
probability and the achieved error (1-fidelity) for fixed Q = 2g, k, = 10g,'jA = 1b = 
0.1(7 and A = 2Qg and a waiting time of T = 100/(7 against the detector efficiency. 
We observe, as expected, a weak linear reduction in the state fidelity see that the 
success probability scales by a factor 77. 

5.4.3 Dark counts 

A further experimental imperfection which must be considered is the presence of 
"dark counts", i.e. when the detector fires although no light is incident upon it. 
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Clearly, this will degrade to some extent the performance of all schemes which rely 
on a single detector click to generate an entangled state, and will lead, in general, 
to a loss in fidelity of the state produced. However, in the present scheme, the 
time-window in which a click due to a photon should occur is far shorter than the 
mean time between dark counts. For example, in [ButtlerOO] a dark count rate of 
approximately 1400 s~^ is reported, thus the mean time between dark counts is on 
the order of ms. In the optical regime, the atom-cavity coupling g, detuning A, 
cavity decay rate k and the coupling with the classical field fl will all be at least on 
the order of MHz [McKeeverOS]. Thus, using equation ()5.24|) one can estimate that 
Tav, the mean time before a proper click occurs in this scheme, is on the order of /is. 
Since the time-window for detection in this scheme can thus be made much smaller 
than the mean time between dark counts, their effect on this scheme can be made 
very small. 

So far we have considered the case where we were faced with the task of entan- 
gling two spatially separated ion-cavity systems. The key ingredient in the detection 
scheme was the beam-splitter that erased the which-path information from the sys- 
tem, so that a photo-detection event could lead to entanglement between the cavities. 
However, the above method could also be used to entangle two ions trapped in a 
single cavity whose decay is monitored by a single photo-detector if the system is 
set-up such that the detection of a photon does not provide any information about 
which ion the photon was emitted from. 

5.5 Conclusion 

In summary, we have presented here a new approach to entangling ions in spa- 
tially separated cavities that, under ideal conditions, allows for the deterministic 
generation of perfect entanglement between individual ions in distant cavities. Em- 
ploying a weak interaction as in the previous chapter provides a robust method of 
state-preparation which is resilient against many experimental imperfections. We 
therefore expect this scheme to work well when operated with current experimental 
setups. The reason why this approach is so effective can be understood as follows. 
Ignoring decoherence and dissipation, the evolution of the system is unitary and can 
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be thus written as an exponential exp[—ixO], where O is a Hermitian operator of 
unit norm. In the regime of weak or brief interactions, x is small, and the operator 
may thus be approximated by 1 — ixO. Then, a measurement takes place which 
projects out one of these terms. With high probability the former term is selected, 
and the system has returned to its previous state, but with a small probability, the 
latter term, which generates the desired entangled state, is realised. By keeping x 
small the actual value of x remains unimportant. While in the first scheme, pre- 
sented in chapter m this led to a low success probability, in the proposal presented 
in this chapter, constant weak driving and continuous detection allow, in the ideal 
scenario, arbitrarily high success probabilities to be achieved. The schemes remain 
robust in the face of decoherence, dissipation and imperfect measurement, as long 
as measurements can still resolve that latter term sufficiently well. Thus, such tech- 
niques could be adapted to many other physical systems, and may be especially 
useful in systems which are hard to control with precision. 
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Part III 

Entanglement Distillation with 
Linear Optics 
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Chapter 6 



Gaussian States and Operations 



In the second half of this thesis, we move from the cavity QED setting to consider 
the quantum states of hght pulses. As we discussed in section 12. (j^ the quantum 
state of a stream of pulses generated by an ideal laser source may be written, 



As we discussed in section 1221 as long as all phase measurements on the light are 
made using reference pulses split from the same source (which is just the approach 
that experimentalists take) the pulses will be indistinguishable from a pure product 
of coherent states, lAe**^). In the following chapters, for mathematical and conceptual 
convenience we will shall use the universal approach taken in quantum optics and 
ignore the uncertainty in the overall phase of the generated light. This is perfectly 
valid as long as a beam of reference pulses split from the source as described in 
section 12.61 is available to the experimenter along with the pulse whose quantum 
state is being manipulated, and we assume the presence of such reference pulses in 
all the following discussion. 

6.1 Gaussian states 

Clearly the quantum state of a light mode can be extremely complicated, as a general 
state of a quantum harmonic oscillator has infinitely many parameters. Fortunately, 




(6.1) 
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there is a class of states which require only a small number of parameters for their 
specification, which also have immense experimental importance, especially in the 
frame-work of linear optics. These are the Gaussian States. Gaussian States are 
states whose characteristic function, introduced in section 12. 4| are Gaussian [Si- 
mon87, Arvind95, Eisert03a] . 

In section 12. 4t we restricted ourselves to single-mode states, but to describe 
entangled states, we need a multi-mode notation. Let us label each mode of the n- 
mode system with a number i, and write the quadrature operators (again introduced 
in section 12. 4|) . for the ith mode Xj and Pj. It is convenient to represent these 
operators in terms of a column vector^ R = {Xi, Pi, ... , P„), and likewise form 
a vector of the multi-mode phase-space variables ^ = (Xi, Pi, . . . , Xn, Pn)- 

The canonical commutation relations can now be expressed 

[Rj,Rk]=i'Lj,k (6.2) 

where E is the symplectic matrix, a direct sum of n anti-symmetric two-by-two 
matrices, 

-§(::)• 

We can thus express the multi-mode Weyl operator as follows 

= e'^^^^ (6.4) 

where we have employed the symplectic product, ^"'"SR between vectors ^ and R. 
We write the characteristic function for an n-mode state p, 

Xp(0=Tr[pW^^] . (6.5) 
The state p may be obtained from its characteristic function as follows 



P = J^ I X,i-OW,d^-^ , (6.6) 



^As with the quadrature operators themselves, conventions for the definition of all of these 
multi-mode constructs, in terms of their normalisation and ordering, vary widely in the literature. 
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which can be verified using the identity Tr[#aW^/3] = (27r)"5^"(Q: + f3) [Cahill69]. 
Gaussian states are all those for which the characteristic function has a Gaussian 
form. This means that the characteristic function may be written^, 

Xp = exph^e^ErE^I + ^d^E^] . (6.7) 

The 2n x 2n matrix F is called the covariance matrix and d is the displacement 
vector. The displacement vector describes the centre of the Gaussian in phase-space, 
and the covariance matrix, naturally, the variances and co-variances. By equating 
partial derivatives of Xpi^f^ + t(^) terms of s and t, in equations ()6.5p and ()6.7|) . 
one obtains expressions for F and d in terms of p. The displacement vector consists 
of nothing more than the expectation values of the quadrature operators, 

dj = Tr[pR,] = {R,) . (6.8) 
The elements of the covariance matrix have the following form, 

F,-fc = Tr[p{(i?, - {R,)), (R, - (i?,))} J 

= 2Re(Tr[p(i?, - d,)){Rk - 4))]) (6-9) 
= 2 Tr [p{Rj - dj){Rk - 4)] ) - i^j,k 

where {A, B}+ = AB — BA is an anti-commutator. We see from this definition 
that F is a real symmetric 2n x 2n matrix. Thus an n-mode state is described 
by 2n + (l/2)(2n)(2n + 1) = 3n + 2n^ parameters. Furthermore, the displacement 
vector contains only quantities local to each mode, which means that all of the 
entanglement properties of the state must only depend on the state's covariance 
matrix. 

It is important to mention that not all possible covariance matrices correspond 
to a physical state. In fact, for p to be Hermitian and normalisable, the covariance 
matrix must satisfy the following inequality 

F + iE>0, (6.10) 
^Again, the normalisations in this are a matter of convention. 
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where the matrix inequahty expression means that the right hand side is positive 
semi-definite. This inequahty also contains the Heisenberg uncertainty relations for 
the quadrature operators. 

Not only do the Gaussian states have a useful theoretical description, they corre- 
spond to states of light of great experimental importance. These include the coher- 
ent states which we introduced in section thermal states of a radiation field and 
the squeezed states [Stoler70, Yuen76], a very important class of non-classical states 
which have been heavily studied in quantum optics (see [Loudon87] or [Teich89] for 
a review), both as single mode and entangled multi-mode states. 

6.1.1 Descriptions of important Gaussian states 

We shall briefiy review how these important states are represented in the covariance 
matrix and displacement vector formalism. The covariance matrix for a single-mode 
coherent state \a) is simply the identity matrix 1 and their displacement vector is 
(Re(\/2a), Ini(A/2a))^. When a radiation mode is in thermal equilibrium with a 
heat bath at temperature T its state is called a Gibbs thermal state, 

p= (l-e-")^e-""|n)(n| (6.11) 

n 

where k = {ksT)^^. Thermal states are Gaussian States with covariance matrix, 
r = coth(/s:/2)l and zero displacement. 

The Heisenberg uncertainty relation for the canonically conjugate quadrature 
operators is AXAP > 1/2. For coherent states, both AX and AP are \/T/2, so 
these are states of minimum uncertainty. To reduce AX further one must increase 
AP or vice versa. Minimum uncertainty states with AX or AP less than a/1/2 are 
called (single-mode) squeezed states. 

It is convenient to express the squeezed states in terms of the single-mode squeez- 
ing operator S{() 

(6.12) 

where C = ^^^'^ is the complex squeezing parameter. Single mode squeezed states 



5(C) = exp 



C"t2 I C -2 

— a' H a 

2 2 
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are therefore expressed 5'(C)|0). The degree of squeezing^ is quantified by r and 
the angle of the quadrature in phase-space which is being squeezed is expressed by 
(f>. The uncertainty of the squeezed quadrature is reduced by a factor e"^ and the 
uncertainty of the anti- squeezed conjugate quadrature increases by a factor e^. This 
is directly captured in the form of the covariance matrix, 

e-^'- 

^ = 7^(0y^■| ^ \n{<f>), (6.13) 



2r 



e 



where 

, , / cos((/)) — sin((/)) , 

^('^^ ^ • ! Z I (^-^^^ 

\ sm(0) cos(0) 

is a phase-space rotation through angle (p. The displacement remains unaffected 
by the squeezing operator and remains zero, although a more general family of 
squeezed states may be constructed by allowing finite displacements as well. Exper- 
imentally, the highest degree of squeezing which can been obtained for a continuous 

laser beam is 7.2dB [Schneider97] corresponding to r = 0.82. The highest degree of 
squeezing recorded for pulsed light is S.ldB or r = 0.36 [Wenger04]. 

So far, 1 have only detailed single-mode Gaussian states. The next case is an 
important example of a two-mode entangled state. It is called a two-mode squeezed 
state and is again best expressed in terms of a two-mode squeezing operator, as 
^2(C)|0)|0), where 

'S'2(C) = ~CcilCl-2 + C*0-lO-2 (6.15) 

and where modes are numbered 1 and 2. We shall express the complex squeezing 

parameter C in terms of its modulus and argument, C = re^'^. For non-zero r this 

state is entangled, and its entanglement increases with increasing r. Its covariance 

^If one interprets the uncertainty in the squeezed quadrature measured in an experiment as 
"noise", then squeezing may be quantified in terms of of dcciBcls. In fact, this is the approach 
taken in most experimental papers. The conversion of a squeezing value d expressed in deciBels to 
the squeezing parameter r proceeds via the following relation r = (ln(10)/20)d w 0.12d. 
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matrix is written, 



r = [7^(0/2) © n{(j)/2)f 



/ 







Sf 





\ 


















Sf 





Cf 
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J 



7e(0/2) © 7^(0/2) (6.16) 



where C,. = cosh(2r) and 5",. = sinh(2r). The displacement vector for this state is 
zero. 

One of the beautiful features of the covariance matrix formalism is the ease with 
which one can identify the reduced state of a sub-system, one merely reads off the 
corresponding sub-matrix. Thus, for example, we see that the reduced state of two 
modes in a two-mode squeezed state is a thermal state with coth(K/2) = cosh(r/2). 
By writing out the state vector of the two-mode squeezed state, we see that it also 
possesses correlations between the photon numbers of the entangled modes. 



^2(re^'^) =sech(r)^[-e*'^tanh(r)]"|n)|n) . 



(6.17) 



n=0 



The "squeezing" in this state is an enhancement in the correlations between the 
quadratures of the two modes. Let us set to zero and consider the uncertainty in 
the sum of the X quadratures for both modes. 



A(Xi + X2) = ^((Xi + X2)2) 



ri,i + r3,3 



+ ri,3 



(6.18) 



We see that, in this case, the sum of the quadratures is squeezed, and similarly one 
can show that the difference of the quadratures is anti-squeezed. When we consider 
the sum and difference of the conjugate quadratures A and P2, we find that the 
opposite occurs, namely the difference is squeezed and the sum is anti-squeezed. 

Two-mode squeezed states may be generated by mixing two single- mode squeezed 
states, with orthogonal squeezed quadratures, at a 50:50 beam-sphtter [Purusawa98, 
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Silberhorn01,Bowen02], or alternatively, generated directly via a non- linear process 
known as parametric down-conversion [Ou92, ZhangOO]. A detailed treatment of 
the theory of squeezed states, both single and multi-mode, can be found in the 
literature [Ma90, Barnett97]. 

6.2 Gaussian operations 

Now that we have examined the most important Gaussian States, we turn our atten- 
tion to the Gaussian operations [Arvind95], namely those which map all Gaussian 
states onto Gaussian States. We shall see that these operations have significant 
experimental importance in the optical setting. We will not discuss the formal de- 
scription of the Gaussian operations in great detail here it is worth briefly describing 
how operations on states translate into operations on covariance matrices and dis- 
placement vectors. We will focus here on unitary operators alone, for an excellent 
review which includes more general operations and measurement, see [Eisert03a]. 

Unitary transformations on Gaussian states correspond to symplectic transfor- 
mations on covariance matrices and displacement vectors. A symplectic transforma- 
tion S is one which preserves the symplectic product between any pairs of vectors, 
and corresponds to a canonical transformation of quadrature operators [Ekert91b], 
preserving the canonical commutation relations. This leads to the condition, 

S^ES = E (6.19) 

where E is the symplectic matrix defined above. Symplectic transformations are rep- 
resented by real matrices with unit determinant. For a given unitary transformation 
on a state U, there is a corresponding unique real symplectic transformation such 
that UW^W = Ws^- It is then straightforward to confirm that when a Gaussian 
state is acted on by U, its displacement vector transforms d i— > Sd and its covariance 
matrix, F i-^ STS^. Once can show (see [Eisert02], [Giedke02] and [Fiurasek02] ) , 
that, in the optical setting, all Gaussian operations can be implemented with beam- 
splitters, phase shifters, squeezers and homodyne measurement. Apart from the 
detection, these are the operations which comprise linear optics. 
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Linear Optics is so-called as it encompasses all operations which can be expressed 
as linear transformations of the creation and annihilation operations. For example, 
a single mode pure state 



|^) = J]a„|n) = J]^(atr|0) (6.20) 

n n 

is transformed by a single mode linear transformation L[-] to 

l^') = E^(^[«']riO)- (6-21) 

These operations may be divided into two categories, passive linear operations, 
which preserve the mean photon number of the state, and active linear operations, 
which do not. 



6.2.1 Passive linear optics 

All passive linear operations can be implemented via networks of phase shifters and 
beam splitters. Phase shifters are single mode operations which transform to 
e^'^d'^ag, where is a particular phase. Thus, for example, the coherent state \a) is 
transformed to the state \e^'^a). Phase shifts are often implemented by passing the 
laser pulse through a crystal with a certain length with a refractive index higher 
than air. The symplectic transformation which the phase shift generates is, as one 
would expect, a rotation in phase space 

^ ^ / cos(0) - sin(0) 
\ sin(0) cos(0) 

The beam-splitter takes two-modes as input and outputs two-modes. It im- 
plements a linear two-mode transformation given here in its most general form 
[Prasad87] 




(6.23) 
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where |Tp + = 1. T (or sometimes |Tp) is referred to as the transmittivity 
and R (or sometimes |-Rp) as the reflectivity. Thus for a 50:50 beam sphtter |Tp = 
\R\^ = 1/2. 

Beam-sphtter operations are typically implemented by allowing the optical pulses 
to meet on a semi-silvered mirror. The phases of T and R may then be chosen by 
appending extra phase shifters to the inputs and/or outputs of the beam splitter, 
as appropriate. In this thesis, we shall assume that T and R are always real, and 
always specify any extra phase shifts. In that case, the symplectic beam splitter 
transformation is 




(6.24) 



where I2 is a two by two identity matrix. 

Approximate displacement operations may be achieved with linear optics by 
mixing the mode to be displaced with an intense coherent state on a beam-splitter 
with near 100% transmittivity. 

6.2.2 Active linear optics 

In contrast to the passive linear operations "squeezing" , the physical realisation of 
the single mode squeezing operator S{(), changes the energy of the mode and thus 
requires driving by an external field. This field b, couples to the mode which is to 
be squeezed a inside a non-linear crystal with the interaction Hamiltonian, 

H = t{xiafb^-X*iaYb) ■ (6.25) 

The driving field is an intense "classical" coherent state and thus, taking the 
classical approximation introduced in section 13. 5t this implements the squeezing 
unitary 

Si2x*(3) = exp [-x*Pia^)' + xP*ia)'] • (6.26) 
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The symplectic transformation for a squeezing operation with real ( — r is, 



A complex squeezing parameter corresponds to extra phase rotations, as above. 
Of course, many approximations have gone into the analysis here, and in physical 
realisations there is both a limit in the degree of squeezing which is attainable, and, 
typically, extra unwanted noise will be added to the state. 

6.2.3 Homodyne measurement 

The statistics of the quadrature operator on a mode may be probed by homodyne de- 
tection [Personick70,Yuen80,Banaszek97]. This is implemented by mixing the mode 
with a phase coherent reference pulse (i.e. from the same source) in an intense co- 
herent state on a beam-splitter, and then measuring the difference in intensity of 
the two output ports. The measurement corresponds to a projection onto the eigen- 
state of X with the measured eigenvalue. By adding phase shifts to the reference 
pulse, other quadratures may be measured. A more elaborate setup using four beam 
splitters, known as eight-port homodyning, leads to projections onto coherent states. 

6.2.4 Non-unitary Gaussian maps 

If one further ingredient is added, the above operations include all possible Gaus- 
sian operations. This is to go beyond merely unitary operations and projections, 
by adding ancilla modes prepared in some fiducial state, allowing joint operations 
between our system and these ancilla modes, and tracing out the ancilla modes af- 
terwards. The general completely positive Gaussian maps which this leads to are 
characterised in [Eisert02, Giedke02, Fiurasek02]. 




(6.27) 
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6.3 Entanglement in Gaussian states 

As one would expect for a system with infinitely many degrees of freedom, in gen- 
eral, quantifying entanglement in infinite- dimensional systems can be difficult. In 
particular, many standard entanglement measures show undesirable characteristics 
such as discontinuities [EisertOSb], and one can show that any state is arbitrarily 
close to a state of infinite entanglement. However, as Eisert and co-authors show, 
if one makes the reasonable demand that all states must be finitely bounded in en- 
ergy, then these problems do not arise, and entanglement measures are much better 
behaved. 

The matter is simpler still for Gaussian states, where, due to their small number 
of parameters, particular entanglement measures can be calculated directly. Addi- 
tionally, since the displacement vector contains purely local properties, the entan- 
glement of a state depends on its covariance matrix alone. 

The question of separability of Gaussian states was first considered independently 
in [DuanOO] and [SimonOO]. In particular, Simon showed that the Positivity of Partial 
Transpose (PPT) criterion is both necessary and sufficient for bi-partite Gaussian 
states. In phase-space, transposition corresponds to time-reversal, momentum is 
reversed, position remains unaltered. Consider a Gaussian state p with covariance 
matrix F in a system of n x m modes. The entanglement between the former n 
and latter m modes may be investigated by considering the partially transposed 
covariance matrix F' where the momenta of the latter m modes have been reversed. 
Thus the PPT criterion becomes the condition that F' fulfills the inequality 

F' + > . (6.28) 

If this condition is not fulfilled, the state p must be entangled. Werner extended 
this result to show than if entanglement between a single mode in a multi-mode 
Gaussian state and the remaining modes are consider this criterion is also necessary 
and sufficient [WernerOl]. 

As we described in section II. 2| one can derive a measure of entanglement from 
this condition, the "logarithmic negativity" , which expresses to how great a degree 
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the PPT inequality is violated. For Gaussian States, this quantity can be simply 
expressed in terms of the symplectic eigenvalues of the covariance matrix. The 
symplectic eigenvalues [Simon94] are the absolute eigenvalues of the matrix EF, 
which always occur in pairs. Any covariance matrix can be brought into a diagonal 
form with these entries by symplectic operations [Williamson36]. For a state p, 
with covariance matrix F, if F' is the partially transposed covariance matrix, the 
logarithmic negativity is defined to be, 

£;Ar = ^Max[-log2(7D,0] (6.29) 

k 

where 7j[. are the symplectic eigenvalues of F' and only one member of each pair of 
eigenvalues is included in the sum. 

Two further important quantities can be expressed concisely in terms of the 
symplectic eigenvalues, the Von Neumann entropy Sy^ — — Tr[plnp] and linear 
entropy Sn^ = 1 - Ti[p'^], 



E 

k 

and 



— log — - — — log 



(6.30) 



Siin = l-l\- = l-\r\-'^ , (6.31) 

where 7^ the symplectic eigenvalues of F and, as above, each eigenvalue pair is 
counted only once. The determinant of F is written |F|. This is equal to the square 
of the product of the symplectic eigenvalues since the determinant is unchanged by 
symplectic (unit determinant) operations. A corollary of this statement is that pure 
states always have covariance matrices with unit determinants. For a full derivation 
of these equations see, for example, [KriigerOl]. The expression for the von Neumann 
entropy allows us, of course, to calculate the entanglement of pure Gaussian states. 
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6.4 Gaussian states in quantum information 

The experimental feasibility of Gaussian operations makes Gaussian states an at- 
tractive setting for implementation of quantum information science protocols. These 
protocols have in the main been proposed for qubit systems, however many may be 
translated to the continuous variable domain. 

Braunstein and Kimble showed how quantum teleportation could be imple- 
mented with linear optics using two-mode squeezed states as an entanglement re- 
source instead of Bell states [Braunstein98b] and their proposal was implemented 
soon afterwards [Furusawa98] . Quantum key distribution (QKD) also translates 
readily to the continuous variable regime. Optical states are particularly suitable for 
this task as they can be easily transported over long distances. QKD schemes based 
on the transmission of coherent states [Grosshans02] , squeezed states [RalphOO] 
and shared two-mode squeezed states [RalphOO] have been proposed and reahsed 
[Grosshans03, Silberhorn02b, Silberhorn02a] . So far only the latter two have been 
proven absolutely secure against all possible eavesdropping attacks [GottesmanOl]. 

Other qubit-based protocols which have been translated into linear optical pro- 
posals on Gaussian states include dense-coding [BraunsteinOO], entanglement swap- 
ping [van LoockOO] and error-correction [Braunstein98a] . An important protocol 
which had not, however, been developed for continuous variable states is entangle- 
ment distillation. In the following chapters we shall describe a new protocol for the 
distillation of entangled Gaussian states. 
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Distillation of Continuous Variable 
Entanglement 



7.1 Distribution of entanglement 

As discussed in chapter [TJ many of the most intriguing features and useful apphca- 
tions of entanglement are encountered when the entangled sub-systems are spatially 
separated. However, entangled states can only be created by global operations on 
the sub-systems, which requires either that the sub-systems interact directly with 
one another, or each interacts directly with some mediating third party, as in the 
protocols described in chapters 0] and El In either case, for entanglement to be estab- 
lished between two distant laboratories, quantum systems must physically traverse 
the distance between them. 

It is unavoidable that during this transmission, some interaction with the "en- 
vironment", which is how we shall describe systems outside of the experimenter's 
control, will occur. This usually leads to the state of the system changing, typically 
becoming more mixed due to entanglement with the environment. This in turn 
degrades the entangled state which is supposed to be being transmitted. 

As an example of this, let us consider a two-mode squeezed state, as introduced 
in section with real squeezing parameter A, generated at a central source of 
pulsed radiation. The two pulses are transmitted separately through optical fibres 
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to two parties, who we shall name Alice and Bob^. 

To a first approximation, an optical fibre acts as an absorbing channel. This is 
a transmission line which absorbs some of the radiation which passes through it. 
A thorough analysis [ScheelOl] has shown that such a channel may be accurately 
modelled in a very simple way. We pretend that the light pulse has been mixed with 
a vacuum mode on a beam splitter of transmittivity T = ^/T. The parameter r 
then corresponds to the transmission coefficient of the channel, expressing the ratio 
between the intensities of transmitted and incident light. Since the beam splitter 
is a Gaussian operation, the absorbing channel preserves the Gaussian nature of 
input states. In the phase-space formalism it is easy to calculate the effect of such a 
channel. The displacement vector is reduced by a factor of i/r, and the covariance 
matrix F transforms to rF + (1 — r)l. 

Let us consider the effect of passing both pulses of a two-mode squeezed state 
with real squeezing parameter r through two channels, each with equal transmis- 
sion coefficients. The covariance matrix is transformed, from the form expressed in 
equation to 

/ l + r(a-l) -rSr 

1 + T{Cr - 1) 

-rSr 1 + T{Cr - 1) 

^0 rSr 1 

where = cosh(2r) and = sinh(2r) as before. We can calculate the logarithmic 
negativity of this state using the method described in the previous section, and find, 

£;^ = -log2[l-r(l-e-20] • (7.2) 

Thus, as expected, while entanglement still increases with increasing r, the absorp- 
tion r < 1 reduces its value. Furthermore, while the entanglement of the pure 

^The use of the names Alice and Bob, instead of the dry labelhng A and B, has been adopted 
from classical communication theory to become ubiquitous in quantum information science. Indeed, 
it is hard to find a paper on quantum key distribution or entanglement distillation where they do 
not appear! 







(7.1) 
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two-mode squeezed state is unbounded, the entanglement of the state after the ab- 
sorbing channels, is bounded by — log2 [1 — r] . 

7.2 Entanglement distillation 

This degradation in the entanglement of the states at their source is a generic feature 
of the decoherence processes which are unavoidable in the transmission of quantum 
states. For the applications of spatially distributed entangled Gaussian states dis- 
cussed above in section 16. H highly entangled states are required. Since, by defini- 
tion, the entanglement of a spatially distributed state cannot be increased by local 
operations and classical communication, another approach is required. This prob- 
lem was first considered for qubit systems, and Bennett and co-workers proposed a 
scheme [Bennett96c] whereby if the two spatially separated parties, Alice and Bob, 
share a supply of many copies of a mixed and weakly entangled state, which had 
a fidelity greater than 1/2 with respect to a maximally entangled state, they could 
distill a small number of copies of states arbitrarily close to a maximally entangled 
state by two-qubit operations, measurement and classical communication. 

Clearly, the quantum systems most appropriate for transmission over long dis- 
tances are light pulses. Light pulses can embody qubits in, for example, their po- 
larisation, and schemes for the implementation of the quantum gates needed for 
entanglement distillation have been proposed [KnillOl]. However, these schemes are 
very technically challenging and in particular, a full scheme would require detectors 
with close to unit efficiency, which are so far unavailable. So far, therefore, only a 
few steps of such a scheme have been implemented [Pan03]. 

Of course, light pulses can also distribute entangled Gaussian states, which we 
have seen in the previous chapter, have many applications in quantum information 
science. It is thus natural to ask whether one can distill Gaussian states. In partic- 
ular, since Gaussian operations are so experimentally accessible, it would be ideal 
to have an entanglement distillation protocol for Gaussian states which employed 
solely these operations. 
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7.3 Distillation with Gaussian operations 

With the formal description of Gaussian operations, which we have outhned above, 
the question of finding a Gaussian entanglement distillation procedure becomes an 
matter of matrix analysis. This was carried out in [Eisert02], [Fiurasek02] and 
[Giedke02]. Their proofs are too technical to sketch here, but, in light of what is 
possible with Gaussian operations alone, the result is surprising. If one acts with 
solely Gaussian operations and Gaussian measurements on a supply of Gaussian 
states it is impossible to distill any greater entanglement than the initial states 
already possess. Any joint operations, measurements and classical communication 
which Alice and Bob implement will leave them with a state with either the same 
or less entanglement than the initial states. 

This seems a disappointing result, but there is a loophole, that means that optical 
distillation of Gaussian states is not entirely ruled out. There is one operation 
available to the linear optics lab, that is not a Gaussian operation. This is photo- 
detection. In the next section, we will show that even very simple protocols using 
passive linear optics and photo-detection can allow Alice and Bob to generate states 
more highly entangled than their initial supply, but the states will no longer be 
Gaussian. 

7.4 Procrustean entanglement enhancement with 
linear optics and photo-detection 

Imagine Alice and Bob share a pure two-mode squeezed state with real squeezing 
parameter r and Bob performs the following operation. He mixes his entangled 
pulse with a single photon on a beam splitter with transmittivity T as illustrated 
in figure 17. and measures with a photo-detector one of the beam splitter outputs. 
Let us assume for now that the photo-detector works with 100% efficiency and can 
resolve any number of photons in the measured mode. It is then a matter of algebra 
to show that after the measurement, if m photons are registered Alice and Bob's 



106 



7.4 Procrustean entanglement enhancement with hnear optics and 

photo-detection 



Alice 



Classical 
Channel 

Source 



Bob 



Photo- 
detector 




11) 



/BearrN 
Splitter 



Photon 
Source 



Figure 7.1: A "Procrustean" linear optical entanglement enhancement protocol. Bob 
mixes his half of a two-mode squeezed state with a single photon on a beam-splitter of 
transmittivity T and places a photon counter at one of its output ports. 

state is collapsed into the following, 

oo 

an\n) A\n + 1 — m) (7.3) 



T.n=m-1 l«nMI 



2 

n=m— 1 



where 



a^(m) = (-tanh(r))"T"-™i?'^~i 

V V"'' / V V"^ ~ / 

(7.4) 

and we have assumed for simplicity that T and R are real, and we define (^) = if 
k > n. 

As this is a pure state, one can immediately calculate the entanglement shared 
by Alice and Bob in terms of the von Neumann entropy of the reduced density 
matrix pA, E^n = — Tr [p a log2 Pa]- The analytic form of this does not give us much 
insight, so we shall plot the entanglement of the state generated when m photons 
are detected, and the probability that this will occur (given a perfect detector) in 
figure 17.21 

This shows that with this very simple protocol employing the non-Gaussian 
operation of photo-detection, a huge gain in entanglement is achievable, if at a very 
low probability. The effect of the beam-splitter and the photo-detector is to single 




m) + T 
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Figure 7.2: The pure state entanglement (a) of the state generated and the event's 
probabihty (b), when a certain photon- number is detected in the simple protocol described 
in section TlAl Alice and Bob share on a initial two-mode squeezed state with squeezing 
parameter tanh(r) =0.1 and a beam-splitter of transmittivity T = tanh(r) =0.1 is used. 
The probabilities which are close to zero and can thus not be read clearly on the graph 
decrease exponentially, e.g. for 2-photons the probability is 3 x 10~^ and for 3-photons the 
probability is 4 x 10~^. For comparison, the pure state entanglement of the input state is 
0.08. 



out a particular decomposition of the input state, of which the part corresponding 
to zero-photons in the detector has a high amplitude and very weakly entanglement, 
and the 1-photon part has low amplitude and high entanglement. For the low values 
of r which are experimentally attainable, the probability of a two-or-more photon 
detection event is low, however, the beam splitter parameters can always be chosen 
such that, as in figure 17^ the one-photon detection event corresponds to the most 
highly entangled state. 

The fact that a degree of entanglement of almost unity is obtained in the above 
example is not accidental, since the value of T was been chosen to ensure that 
the state generated in close to a maximally entangled two qubit state in the |0) 
and |1) photon basis. If one chooses T ^ tanh(r)/A for some value A, then the 
state generated conditional on one photon being detected is approximately |0)|0) — 
A|l)|l), the approximation being more and more valid as the original squeezing r 
decreases. Thus the above method generates approximate Bell states, conditionally 
on a particular measurement outcome being obtained. In the next section we shall 
see an important and relevant use for such states. 

So far we have assumed perfect conditions, which we shall now relax. Firstly, we 
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have assumed that the detector is both 100% efficient and resolves different photon 
numbers. However, since the probabihty of a two-or-more photon detection event 
occurring is so small, when the detector does click it can be assumed that only one 
photon was present, thus the use of an inefficient detector which merely registers 
the presence or absence of photons will only minimally affect the state produced, 
leading to a tiny admixture of the states from these multi-photon detection events. 

Secondly, if Alice also possesses such a photo-detector. Bob does not need a 
single-photon source, which are a great technological challenge. If the pair is willing 
to sacrifice some entangled pairs. Bob need only wait till Alice detects a single 
photon on her side of a pair and this will project Bob's counterpart into a single 
photon which he can use. 

Most importantly, if the protocol were realised between distant laboratories, the 
light pulses will have passed to Alice and Bob via a channel which will at the very 
least absorb some of the signal decreasing the entanglement and increasing the en- 
tropy of the initial state available to Alice and Bob. Thus it is natural to ask whether 
this protocol can give a similar boost in entanglement for a mixed input state. We 
consider here the state described by the covariance matrix in equation ()7.1|) . The 
matrix elements of the state generated by the Procrustean protocol, in the case 
when one and only one photon is incident on the detector may be calculated analyt- 
ically, and the entanglement of the state, in terms of its logarithmic negativity thus 
calculated. In figure 177^ we plot the entanglement of the state before and after a suc- 
cessful implementation of the protocol. We see that even for these mixed states, this 
extremely simple procedure can provide a significant increase in the entanglement. 

It will be convenient for the next chapter to write out the zero-photon and one- 
photon matrix elements of the (unnormalised) state p generated here. 



Po,o,o,o 



1 



(7.5) 



can be set for convenience since Po,o,o,o > 0. 



Pi, 1,0,0 — Po,o,i,i — '^Po,o,o,o 



(7.6) 
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Figure 7.3: In this figure we see the result of the procrustean entanglement concen- 
tration protocol on a mixed input state which is the result of a two-mode squeezed state 
with squeezing tanh(r) = 0.1 being passed through absorbing channels with transmission 
coefficient r. The beam splitter transmittivity is chosen to be T = 0.1, the same as the 
plot in figure [721 In plot (a) the entanglement of the generated state is plotted alongside 
the entanglement of the input state (dotted line) as a function of r. In plot (b) we see 
that the success probability has decreases weakly and linearly with decreasing r. 

Pi,i,i,i = AVoao.o (1 + - 1) tanh2(r)) (7.7) 

Po,i,o,i = eAVo,o,o,o (l - - 1) tanh^(r)) (7.8) 

where 

(2T^-l)tanh(r)r 

r(tanh2(r)(r- 1)2) - 1 ^ ' ' 

and e = (1 — r)/r. All other zero and one photon matrix elements are zero. Note 
that when r is very small, these matrix elements depend only on the two parameters 
e and A, and furthermore for any given values of r and r, A may be set by choosing 
the appropriate value of T. 

The first investigation of local measurement-based protocols for probabilistically 
increasing the entanglement of a state can be found in [Bennett96b] who coined 
the term "Procrustean" for simple measurement-based protocols on a single copy 
of the state^. The protocols described here were introduced by us in [BrowneOSa] 
^Procrustes was a villain in Greek legend, who would invite travellers to sleep in his iron bed, 



110 



7.4 Procrustean entanglement enhancement with hnear optics and 

photo-detection 



and [Eisert04], but a more comprehensive investigation of similar ideas can be found 
in [ThewOla, ThewOlb]. Since any kind of entanglement concentration or distilla- 
tion is impossible with Gaussian operations, it is perhaps surprising how well these 
very simple Procrustean protocols with a single photo-detection as their only non- 
Gaussian element can enhance the entanglement of their input states, both pure 
and mixed. Theses states however, are manifestly non-Gaussian and not useful for 
the applications described in section 16 .41 It would therefore appear, at first sight, 
that these protocols still do not aid us in the search for an entanglement distillation 
procedure for Gaussian states. In the next chapter, however, we introduce a simple 
linear optical protocol called "Gaussification" . We will show that this allows Alice 
and Bob, via local linear optics operations, measurements and classical communica- 
tion only, to convert entangled non-Gaussian states into Gaussian states and also, 
often, further increase the entanglement in the process. 



promising that it would be a perfect size for them. However, this perfect fit was achieved by shorter 
guests being stretched or chopping off the legs of taller ones! The entanglement concentration 
methods here are Procrustean, in the sense that the projection, corresponding to the successful 
measurement outcome, "chops off" the more weakly entangled portion of the state leaving a more 
highly entangled state behind. 



Ill 



Chapter 8 

"Gaussification" with Linear 
Optics 

8.1 Introduction 

In this section we introduce an iterative protocol, that we shall call "Gaussifica- 
tion" [BrowneOSa, Eisert04] . This protocol acts on a supply of many copies of a 
non-Gaussian two-mode entangled state, under conditions which we describe below. 
After several non-deterministic iterations, it produces a small number of states that 
are arbitrarily close to Gaussian states and which are often more entangled than the 
input. First we will outline the protocol and then, under the assumption that all 
operations are ideal, we will consider analytically all its important properties. In 
the final section of this chapter we will model how the procedure would be affected 
by the most important experimental imperfections of a real implementation of the 
scheme. 

Each step of the protocol proceeds as illustrated in figure 18.11 It is an iterative 
procedure. Each iteration consists of 

1. In the first instance, Alice and Bob act on shared two-mode entangled states 
supplied from a source. In later iterations they use pairs of states generated 
by a previous successful iteration. 

2. They each mix their half of two copies of the state on a 50:50 beam splitter. 
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Channel 
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Figure 8.1: A single iteration of the "Gaussification" protocol. 

3. They measure one of the outputs with a detector which distinguishes the 
presence from the absence of photons. 

4. If no photons are registered, the iteration was a success, and the resultant 
state is saved for the next iteration. 

Before we investigate the properties of this protocol for arbitrary input states, 
let us illustrate it using the simple input state \'ip^'^^) = |00) + ?7|11). We label the 
two entangled modes 1 and 2 and use the labels A and B to denote Alice and Bob. 
Thus the initial state of the two copies of \ip^^^) is 



|V>(0)) ® |7/>(o)) =|0)^i|0)bi|0)a2|0)b2 + ^|1)ai|1)bi|0)a2|0)b2 

+ v\0)ai\0)bi\1)a2\1)b2 + V^\'^)ai\1)bi\1)a2\1)b2 . 

The 50:50 beam splitter transformation transforms the state. 



(8.1) 



■^\0)ai\0)a2\0)bi\0)b2 + |(|1)ai|0)a2 + |0)ai|1)a2) (|l)m|0)B2 + |0)bi|1)b2) 
+ |(-|1)ai|0)a2 + |0)ai|1)a2) (-|1)bi|0)b2 + |0)m|l)B2) 



+ ^(-|2)^i|0)^2 + |0)ai|2)a2) (-|2)m|0)B2 + |0)bi|2)s2) 



(8.2) 



so that, conditional on a measurement of zero photons in modes A2 and B2, the 
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Figure 8.2: This figure depicts the logarithmic negativity Ej\/ as a function of r] for, i) the 
state \ip) = + ??^)(|0, 0) +r]\l, 1)) used as input supply for the Gaussification (dotted 

line), ii) the state generated by a single (dashed) and two iterations (wider dashing) of 
the protocol and iii) the limiting two-mode squeezed state which, in the ideal case, will 
be reached after many iterations. We use the logarithmic negativity here for ease of 
comparison with the mixed state cases later in the chapter. 

state is transformed to 

2 

= |0)ai|0)bi + v\'^)ai\1)bi + j\2)ai\2)bi • (8.3) 

The probability of obtaining these measurement outcomes is (1 + r/^ + ?7^/4)/(l + 
?7^)^, which for rj = 0.5 is approximately 0.8. We observe that the state remains a 
superposition of states of equal photon number and also that the ratio between the 
coefficients of |00) and |11) has remained constant. It is simple to calculate the state 
produced after a further iteration. In figures IH^ and IH7^ we plot the entanglement of 
the state produced after the 1st and 2nd iterations, and the success probability each 
time. We see that for all values of t] each iteration generates a state with a higher 
degree of entanglement. Furthermore, the success probability for each iteration for 
all values of r] is high, staying between 25% and 100%. 

To see if these features apply to other input states, we now consider inputs of the 
more general form |'?/'(°)) = ^^nCtn^ln^n). States generated by successful iterations 
retain this form. Let us introduce the notation l^/i^*)) = Xln '^"''l'^' state 
generated after i successful iterations of the Gaussification protocol. By writing out 
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Figure 8.3: Here we see the success probability P for the first (dashed) and second (wider 
dashing) iterations of the protocol with input supply = y^l/(l + r/^)(|0, 0) + r]\l, 1)). 

the action of the beam spHtter and measurement on this general state one finds a 
recurrence relation relating the coefficients of the {i + l)th states to the ith, 



a.. 



(i+l) 



r=0 



r ' " 



U) 



This relation allows us to investigate the effect of many iterations of the protocol. 
First, we can ask which states remain unchanged by the scheme. Let us label such 
a state |'?/')*^°°) and its coefficients a{f^ . These must be solutions to the equations 



2-"E 



r=0 



n 



a 



^.5) 



for all n. We solve these equations inductively. The 0th equation gives us ag°°^ 



[a 



(oo)\2 



, therefore a, 



(oo) 



1. Of course, since the states are not normalised, this 



alone does not specify the state at all. Setting 



(oo) 



1, the n = 1 equation 



(oo) 

gives us a} 



a^i°\ We therefore have a free choice for this parameter, and 
set it to Proceeding we find that aj"""* = and a^^^ = and this pattern 
continues. In fact the general solution to equations ()8.5p is = C", which can be 
verified immediately by substitution. This state is nothing other than the two-mode 
squeezed state with squeezing tanh(r) = 
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Since this family of states are the only fixed points of the iteration, starting with a 
supply of states in the form, lip^^^) = an^\n, n) will lead after many iterations to 
a two-mode squeezed state, a Gaussian state. We can investigate the convergence via 
recurrence relations ()8.4j) . For convenience, let us adopt a normalisation convention 
where a^^^ = 1. As we have already seen, this means that Oq'* remains at this 
value. Similarly, in this normalisation convention a^*^ also remains at its initial 
value. We prove the convergence of the remaining coefficients by induction. We 
first postulate that each coefficient will converge to the value a^n^\ and for that aim 
introduce deviations form this final value 5n^ = an^ — ai"""*. We can now rewrite 
equation ()8.4|) 




Let us take the limit that i goes to infinity and assume that all coefficients ain for 
m < n converge to a^"* in this limit. 

lim = 2(^-") lim (6li^) (8.7) 

Thus as long as 2^~" < 1, that is as long as n > 1, if all coefficients am converge 
to for m < n, am will converge to a^^ in the hmit that i tends to infinity. 
Since we have already shown that a^^ and ai^ converge, the convergence of the 
other coefficients follows by induction. This means that formal convergence in the 
coefficients occurs for all input states. However not all these final states correspond 
to physical states. In particular, when ( > 1, the state Xln C"'!'^? '"') is unnormal- 
isable, and therefore not a physical state. Thus we conclude the following; for an 
large input supply of states in the form |-?/;(°)) = X^n '^"^''l'^; repeated iteration of 
successful runs of the Gaussification protocol lead to the generation of a small num- 
ber of two-mode states with squeezing tanh(r) = —a^^^/a^^ as long as q;[*V'^o'' < ^i 
otherwise the procedure does not converge. 

We can illustrate the way the state converges toward Gaussian form by plotting 
the Wigner function of the single-mode states of the reduced state of one mode alone 
for the initial state with parameter t] set to 0.6, and the state after one and two 



116 



8.2 General input states 



steps. The Wigner function is the Fourier transform of the characteristic function, 



Figure IH3] shows the Wigner function of the Bob's reduced state initially and after 
one and two iterations of the protocol. While initially, the Wigner function is far 
from being a Gaussian in phase space, its non-Gaussian features disappear quickly. 

8.2 General input states 

In the previous section we illustrated the Gaussification procedure with simple pure 
states as input. Now we turn to input states which are completely general pure or 
mixed states p^^-* . We will proceed in the same line of analysis which was applied in 
the previous section. We will thus prove that; 

1. All two- mode Gaussian states with zero displacement vector are invariant un- 
der the Gaussification iteration. 

2. Beginning with an arbitrary input supply, the protocol will either lead to 
convergence to a particular Gaussian state or it will not converge. 

3. There are analytic forms for the convergence and the limiting states for a given 
input p^^^ in terms of its matrix elements. 

The proof that all two-mode Gaussian states are invariant can be achieved in a 
straight-forward manner, using the covariance matrix formalism. Let the covariance 
matrix of p^'^^ be F. The covariance matrix of Alice and Bob's pair of two-mode 
states will be F © F, where the modes are ordered Al, Bl, A2, B2. The symplectic 
50:50 beam splitter operation S is 






(8.9) 
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Figure 8.4: The Wigncr function of the reduced state of one system (a) initially and 
after (b) one and (c) two steps for the initial pure state \ip) oc |0, 0) + 0.6|1, 1). 

Since the beam splitter operation will be applied to modes in pairs Al, A2 and 
Bl, B2 (not the above order) we also need a swap operation A 




(8.10) 
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Thus the apphcation of beam sphtters to can be expressed 

r e r ^ (I2 e A e i2)(-5 e -5)(i2 e A e i2)r e r(i2 e a e i2)(,5 e sf{i2 e a e I2) 
= rer . 

(8.11) 

Thus the action of simultaneously mixing modes Al and A2, and Bl and B2 leaves 
the covariance matrix of the states invariant. 

Now let us look at how displacement vector is transformed. Let d be the dis- 
placement vector of p^^^ . The beam splitters transform d thus, 

d e d ^ (I2 © A © l2)(5 © 5)(l2 © A © l2)d © d = © V2d (8.12) 

where is the zero vector. After the measurement of modes A2 and B2, the re- 
sultant state is a Gaussian state with the same covariance matrix as and zero 
displacement vector. This proves the first of our assertions above. 

To analyse the convergence properties of the protocol for a general input we 
proceed as above, constructing a recurrence relation for the matrix elements of the 
state p^^^ produced after i successful iterations; 

abed 

^ V W V M''*'''''"^^*^ n^^ (H ISl 

t^a,b;c,d / j / j / j / j -^"a,b:c.d t^s,t;n,mh'a-s,b-t;c-n,d-m ' yj-^uj 
s=0 t=0 n=0 m=0 

with coefficients M^'^.'"'^ given by 

jijs,t;n,m _ r)-(a+b+c+d)/2 f _t \{a+b+c+d)—{s+t+n+m) 
^^^a,b;c,d ~ ^ \ ^) 

We first observe that all matrix elements p^*), ^ ^ where (a -|- 6 -|- c + d) is odd will 
vanish after the first iteration, since the terms in the sum of equal amplitude pick up 
opposite signs and cancel out. Following the notation which we used above, we label 
a state which is invariant under this map We find that either PoAO,o — ^ 

Po'^Q Q = 0. In this latter case, all other matrix elements are also zero, so we ignore 
this trivial non physical case. If we normalise the state such that PoAO,o ~ 
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following six coefficients of p^°°^ are completely free; Pi°^iO) Po°i oi' Pi'^oi' P2^oo' 
PiXOfi- One can then show by induction that once these 10 parameters (since two of 
the elements are real) are chosen, all other matrix elements are already determined. 
In fact Pa,b,c,d is a polynomial in these elements of degree D < Max{a, b, c, d}. 

We now need to show that all of these fixed points correspond, if they correspond 
to physical states, to zero-displacement Gaussian states whose covariance matrices 
are uniquely specified by these 10 parameters. To do this, we calculate the relevant 
matrix elements for a general zero-displacement Gaussian state. We start with 
equation ()6.6|) . substituting the relevant characteristic function, 

p = / eM~\e^ri:^m^d'-^ ■ (8.15) 

The matrix elements of the Weyl operator can be calculated via its normal ordered 
form [Cahill69]. This gives the following, 

{m\W{^)\n) = 

where are associated Laguerre polynomials. 

We can thus write the matrix element Pa,b,c,d = {(^, b\p\c, d) 



+ ^2 



Pa,b,c,d 



(2vr) 



c\d\ 
d6! 



c?^^ exp 



-^^^s(r + i)s^^ 



6 + ^6 

V2 





-d 

jia-c) 




T ib-d) 
^d 






c 


2 




2 



(8.17) 



We ffist note that all matrix elements Pa,b,c,d, where (a — c) + (a — d) is odd (or 
equivalently {a + b + c + d) is odd) must vanish as the integrand is then an odd 
function. The element Po,o,o,o may be found by Gaussian integration. 



Po,o,o,o 



d^C, exp 



-^^^s(r + i)s^^ 
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The integral converges if F + 1 is positive, and by recognising that the integral 
is invariant under symplectic transformations of F + 1, in particular symplectic 
diagonalisation, one can reduce it to a product of standard Gaussian integrals to 
give, 

4 4 

Po,o,o,o = = — — —T (8-19) 

7i7i |F+1|2 

where 71 and 72 are symplectic eigenvalues of F + 1 and |r + 1| is its determinant. 
We have used here the fact that the determinant of a matrix is invariant under 
symplectic operations. To calculate the remaining matrix elements, we multiply 
out the multinomials in the expression and integrate each individually. Using the 
identity 



-ax +jx 



dj 



(8.20) 



3=0 



each term can be re-written as a differential. We use the following identity to 
"complete the square" 

-^^^E(F + = -^^'^S(F + 1)E^£' + j^E(F + l)-iE^j (8.21) 

where = (£ — 2E(F+l)~^E^j) and thus each multinomial term takes the following 
form. 



^ J Aexp[-l^^E(F+l)E^^]ef^e2'^e3'^e^ 

" P)ki ak2 p)k3 ak4 1 r 1 

yAexp[-i^^E(F + l)E^^+fe] 



df,' df^ dj^^ djt 



(8.22) 



|F+ 11 



Qki Qk2 Qk-i Qki 



exp[j^E(F+l)-^E^j 



Each matrix element can now be calculated by expanding its multinomials and 
differentiating as appropriate. Introducing for convenience the matrix B — (F+l)~^, 
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we find that 



0^1,0,1,0 



1 — — -82,2, 



O"0,l,0,l 



1 ~ s — -84,4, 



1,0,0 



C"l,0,0,l 



0^2,0,0,0 



C"0,2,0,0 



— -Bl,3 — -82,4 + «(-Bi,4 — -82,3), 
2-i/2(-5i,i + 52,2-2^5i,2), 

2-l/2(_53,3 + 54,4 -2^53,4), 
— Bi^3 + i?2,4 — "^(-81,4 + -B2,3)- 



(8.23) 



where aa,b,c,d = Pa,fe,c,d/po,o,o,o- 

Our analysis is now almost complete. We showed above that the most general 
states which are invariant under an iteration of the Gaussification protocol have 
null matrix elements Pa,b,c,d when {a + b + c + d) and are completely specified by 
ten free parameters, corresponding to the six matrix elements highlighted. We now 
see that for any choice of these parameters there exists a "state" with a Gaussian 
covariance matrix F = B^^ — 1 (although it may not be a physical state satisfying 
the condition F > — iS), as long as B^^ exists and is positive. This is the case as 
long as \B\ 7^ 0. However, if we recall that po,o,o,o = 4|i?|^/^, we see that |-B| is only 
zero when po,o,o.o is zero, which is only the case when the unphysical null 

state, so we can discount this. Thus we have shown that the only states invariant 
under an iteration of the protocol are Gaussian states with zero displacement. 

We now consider the criteria for convergence. First we ask, given a particular 
starting state to what state will it converge? Labelling the matrix elements of the 
starting state p^^l ^ ^ and choosing a normalisation such that p^l q g = 1, we find from 
equation ()8.4j) that the six matrix elements highlighted above converge to their final 
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value after a single iteration; 



(i) 
Pl,0,l,0 


— Pi, 0,1,0 


~ Pl,0,l,0 


IPi,o,o,ol 


(i) 
Po, 1,0,1 


~ Po, 1,0,1 


(0) 

~ Po, 1,0,1 


IPo,i,o,ol 


{i) 
Pl,0,0,l 


~ Pl,0,0,l 


~ Pl,0,0,l 


(0) (0) 

Pi,o,o,oPo,o,o,i 


(0 
P2,0,0,0 


~ P2,0,0,0 


(0) 

~ P2,0,0,0 




(i) 
Po,2,0,0 


~ P0,2,0,0 


(0) 

~ Po,2,0,0 


-(1/V2)(pS,o,o)' 


(*) 
Pi, 1,0,0 


~ Pi, 1,0,0 


~ Pi, 1,0,0 


(0) (0) 

Pi,o,o,oPo,i,o,o • 



One can show by induction in an analogous way as above for the simple pure 
state case that the states converge matrix-elementwise after many iterations to the 
Gaussian state with covariance matrix — 1, where B is calculated from the 
matrix elements of p^^-* via the inverse of the relations in equations ()8.23|) . Although 
the matrix elements always converge to their respective final values, the state as a 
whole can only converge if B^^ — 1 corresponds to a physical covariance matrix, 
that is B-^ - 1 > -iS. 

Thus we have proved all the assertions made at the start of this section. Namely 
that the only states invariant under the protocol are Gaussian states centred in phase 
space and that for any given starting supply p^^\ whether or not, and to which state 
it will converge can be calculated analytically from just a small number of matrix 
elements of p''°^ 

8.3 Examples 

8.3.1 An example mixed input 

In chapter 13 we considered the states generated by a simple optical Procrustean en- 
tanglement concentration protocol. In particular we considered the states generated 
when the input state was a two-mode squeezed state where both modes had passed 
through absorbing channels with the same transmission coefficient r. In the limit 
that the initial squeezing was small, this process generates a family of entangled 
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mixed states described by two real parameters A, which can be take any real value 
between -1 and 1, by setting the beam splitter transmittivity in the Procrustean 
protocol to an appropriate value, and r with the following matrix elements 

_ 1 

- l + P + A2(l-r)/r ^^-^^^ 
Pi, 1,0,0 = Po,o,i,i = '^Po,o,o,o (8.26) 



Pi, 1,1,1 — A po,o,o,o (8.27) 
1 — r 2 

Po,i.o,i = A po,o,o.o • (8.28) 

T 

Let us investigate how such a state is transformed by iterations of the Gaussifica- 
tion protocol. Firstly, by applying the procedure describe in section IH^ we calculate 
the covariance matrix F of the Gaussian state to which the state element-wise con- 
verges 



r- A2 



r + A2(2r-l) 

r + A2(2r-r 
2Ar 
-2Ar 



2Ar 






\ 

-2Ar 


A2 + r y 



(8.29) 



This covariance matrix only corresponds to a physical state fulfilling F + zS > if 
the condition r > A^ is met, thus state- wise convergence can only occur when this is 
fulfilled. In figure IH31 we plot, as a function of A, the logarithmic negativity, r = 0.5 
of the input state, the state generated by one and two iterations of the protocol and 
the limiting state to which the protocol is converging. We see a marked increase 
in the entanglement after each iteration, in spite of the fact that the input states 
are not pure. The entropy of these states is depicted in figure IH^ We see that the 
entropy of the final state is, in fact, higher than the starting state, and the entropy 
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Figure 8.5: This figure depicts the logarithmic negativity Ej^ as a function of A for, 
i) the initial state p specified in equations (|8.25j) to (|8.28|) where we set r = 0.5, ii) the 
state generated by a single (dashed) and two iterations (wider dashing) of the protocol 
and iii) the limiting state (solid line). Note that the limiting value only exists for values 
of A < -y/r, since above this value the protocol does not converge. 

increases with each iteration accordingly. In figure IHTfl we see the success probability 
for each iteration. In particular this probability is greater than 50% for most values 
of A below the convergence threshold. Comparing this plot to figure 18.31 indicates 
that the success probability is little affected by the mixedness of the input state. 

We see a simultaneous increase in both the entanglement and the von Neumann 
entropy. This is in contrast to entanglement distillation protocols for qubits, where 
an increase in entanglement is usually accompanied by a decrease in entropy, since for 
any given degree of entropy of the state there is a maximum amount of entanglement 
which the system may possess [Verstraete01,Munro01]. This is also true for Gaussian 
states [Adesso04], although not for infinite dimensional states in general, but for 
a given entropy the maximally allowed entanglement is much larger [Adesso04]. 
This explains the seemingly counter-intuitive behaviour seen here. Of course, for 
application of entangled states to quantum key distribution, for the highest security 
and key rate, pure entangled states are desirable. In the next example we investigate 
the criteria for the convergence of the protocol towards a pure state, and give an 
example to illustrate this. 
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Figure 8.6: Here we plot the von Neumann entropy S'vN for the initial (dotted), first 
(dashed) and second (wider dashing) iterations and the limiting state (solid line) for the 
input state p specified in equations ()8.25|) to ()8.28|) . We see that the entropy of the limiting 
state is non-zero and this is reflected in entropy increases at each iteration. 
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Figure 8.7: Here we see the success probability P for the first (dashed) and second 
(wider dashing) iterations of the protocol with input supply p specified in equations H8.25() 

to ^2B- 

8.3.2 Convergence to a pure state 

It is a straightforward consequence of the results in section 18.21 that certain mixed 
input state supplies will converge to a pure Gaussian state. Any pure two-mode 
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centred Gaussian state p may be written OS'i(Ci)5'2(C2)|0, 0) where O is a general 
passive linear two-mode transformation and Si{Q) is a single mode squeezing op- 
eration [Wolf03]. Since the squeezing operators are quadratic in annihilation and 
creation operators, matrix elements pi,o,i,o, Po,i,o,i and pi,o,o,i are zero. Using the 
map between these matrix elements and the covariance matrix F introduced in sec- 
tion |Hi21 we can verify that these conditions alone suffice to ensure that F has unit 
determinant. Thus these matrix elements being zero is a necessary and sufficient 
condition ensure that p is pure. 

Equations ()8.24|) and this condition leads directly to the following conditions on 
the input supply for convergence to a pure Gaussian state, 



.(0) _ I (0) ,2/ (0) 

A^i,o,i,o ~ Ia^i,o,o,oI /a^o,o,o,o 

pS,o,i = IpSaoIVp&.o (8-30) 

(0) _ (0) (0) , (0) 

Pi,o,o,i ~ Pi,o,o,oPo,o,o,i/Po,o,o,o 

provided po,o,o,o and the conditions of convergence described in section 18.21 are ful- 
filled. 

Let us consider an simple example of a mixed input state which fulfills these 
criteria. 

Po,o,o,o = 1/(1 + Pi,i,o,o = Po,o,i,i = ^/(2 + 2£^) 

Pi,i,i,i = ^V(l + ^') (8.31) 

with e G [0, 1), and Pa,b,c,d = otherwise. 

The limiting state, to which repeated iterations of the protocol on this input will 
tend, is a pure two-mode squeezed state with squeezing tanh(r) = —e/2. How the 
degrees of entanglement and entropy evolve under repeated iteration of the protocol 
on this state are plotted in figures lH^ and lH^ Although for more highly mixed states 
the entropy increases after the first iteration, we see that after repeated iterations 
the entropy decreases for all values of e, as it tends to zero. 

Thus we have demonstrated that under ideal experimental conditions this Gaus- 
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Figure 8.8: This figure plots the logarithmic negativity Ej^ as a function of e for, i) the 
initial state p specified in equation (|8.:-{H) ii) the state generated by a single (dashed) and 
two (wider dashing) iterations iii) the limiting pure two-mode squeezed state (solid line). 
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Figure 8.9: Here we plot the von Neumann entropy 5vN for the initial (dotted), first 
(dashed) and second (wider dashing) and fourth (widest dashing) iterations of the protocol 
on state p described in equation (|8.,31|) . The limiting state is a pure two-mode squeezed 
state, the entropy of which vanished. 

sification protocol to give a scheme which can distill Gaussian states of higher entan- 
glement, and for certain input states, of lower entropy from a non-Gaussian supply. 
Combined with a Procrustean entanglement concentration step, as discussed in sec- 
tion E31 this gives an entanglement distillation procedure for Gaussian states based 
on realistic operations, the only one proposed to date for the continuous variable 



128 



8.3 Examples 

regime. In the next section we shall continue our analysis of the Gaussification 
protocol by considering how well it might perform under the imperfect operations 
which would be implemented in a real laboratory. 
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8.4 Imperfect devices 

In this section we show how the most hkely experimental imperfections affect the 
performance of the Gaussification protocol, focusing on its ability to increase entan- 
glement. The importance of various imperfections will depend upon the complexity 
of the implementation. Here we focus on two regimes; firstly the minimal setup im- 
plementing a small number of iterations in one shot, and secondly a full scheme of 
many iterations which would require the feeding forward of measurement outcomes 
and the storage of pulses in between iteration rounds. In the minimal setup, there 
is one obvious imperfection which cannot with current technology be surmounted. 
That is the fact that photo detectors have finite efficiency, and this is the first 
imperfection which we will consider. 

For a full implementation, storing the light pulses in between iterations will 
clearly cause some degradation in their states, and it is these effects which we will 
focus on in section 18.4.21 In section 18.4.71 we shall discuss other challenges which 
could affect the scalability of the scheme, such as the difficulty of mode-matching 
systems with many optical elements. 

8.4.1 Inefficient detection 

An important experimental imperfection for all implementations is that no photon 
detector exists which functions with unit efficiency. With current technology the 
best detection efficiencies at the frequency ranges typically employed in experiment 
are between 50% and 60%. This means that for every photon incident upon the 
detector there is a finite probability rj that its presence will be registered. Current 
photo-detectors are also unable to distinguish between different numbers of photons 
incident in a single pulse, but this is not a problem for the Gaussification protocol 
which only requires that vacuum is distinguished from the presence of one or more 
photons. However, with today's inefficient detectors even this ability is clearly di- 
minished and one would expect that a protocol which relies on detecting the vacuum 
outcome would perform very badly if inefficient detectors are employed. 

We model an imperfect detector as follows, the usual projection operator onto 
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the vacuum is replaced by 



\o){o\ ^Y.(^-vm{k\, 



.32) 



k=0 



where rj e [0, 1] is the detector efficiency. This formula is equivalent to a beam- 
splitter of transmittivity T = ^ placed in front of a detector with unit efficiency. 

We investigate how the protocol with imperfect detectors by constructing, as 
above, a recurrence relation for the matrix elements of the state p*^*-* after the ith 
repetition of the protocol. We find, after some rather tortuous algebraic manipula- 
tion, that 



where 



oo A+k B+l C+k D+l 

_ 

,c,D - 2^ 2^ 

k,l=0 a=0 6=0 c=0 d=0 



X 



Pa,b 

min{a,j4} min{6,i3} min{c,C} min{d,D} 

E E E E 

s=0 u=0 s'=0 u'=0 



jyra,b,c,d,s,u,s',u' (i) /-x 
^ ^^A,B,C,D h'a,b,c,d J A+k-a,B+l-b,C+k-c,D+l-d 



N' 



a,b,c,d,s,u,s' ,u' 



2-{A+B+C+D+2k+2l)/2 



X 



X 



X 



X 



X 



',A-a+B-b+C-c+D-d+2k+2l 



k+l 



a\(A\lA + k-a\f k 
sj\sj\ A - s J\a-s 

b\/B\/B + l-b\/ I 
uj \ uj \ B — u J \b — u 

c\ fC\ fC + k- c\ f k 



s' 



C~s' 



c - s' 



u'J [u'J [ D-u' )\d-u' 



1/2 



1/2 



.1/2 



.1/2 



(8.33) 



^.34) 



The case rj = 1 describes a detector with unit efficiency, and this expression reduces 
to the one described in the previous chapter. 

To illustrate the performance of the protocol under various conditions we shall 
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use two particular input states as model inputs throughout this chapter. The first is 
the pure state \ipa) = (1/Vl-25)(|0, 0) + 0.5|1, 1)) and the second is the mixed state 
p, = (1/1.5)(|0,0)(0,0| + 0.5(|0,0)(1,1| + |1,1)(0,0|) + 0.25(|1,1)(1,1| + |0,1)(0,1|). 
This is the state in equations ()8.25|1 to ()8.28j) with parameters A = 0.5 and r = 
0.5. Both of these states can be generated by the Procrustean protocol described 
in section 17.41 from a supply of two-mode squeezed states, for \ipa) having been 
distributed through two non-absorbing channel, and for through two channels 
which absorb 50% incident intensity. 

In figures IHTTUl and IS. 1 1 1 we see logarithmic negativity of the input state, \ipa) and 
Pb respectively, and the states generated by one and two iterations of the protocol. 
Remarkably, the performance of the protocol is not greatly degraded by the use of 
imperfect detectors. The only effect we see is a gentle decline in the entanglement 
with decreasing efficiency. This decrease is more severe for p^, and strongest for its 
second iteration. This may be attributed to the significant increase in entropy of 
the state at the first step. In fact, there is a trade-off here between the detection 
efficiency and the mixedness of the input upon which the protocol can still increase 
entanglement. This is seen clearly in figure 18.121 where we plot the logarithmic 
negativity for input states and first iteration, for the example mixed input state p 
defined in equations ()8.25j) to ()8.28|1 with parameters A = 0.7 and varying r. The 
mixedness of the input state increases with decreasing r and we see correspondingly 
that a higher detector efficiency is required to achieve any gain in entanglement. 

Let us discuss why the scheme still functions well with imperfect detectors even 
though it relies on vacuum detections. One important reason is that in the input 
states we have considered, the mean photon number of the pairs of pulses is low, so 
when one detector does fail, it is not very likely that a photon is present to be missed 
by the detector. This feature will be generic, as the scheme is designed for, and only 
converges for, input states which are mostly vacuum. Other symmetries in the input 
states we have chosen enhance this effect. For example, with the input state \ipa), 
there are correlations in the photon numbers of the detected modes which mean that 
errors caused by only one detector failing are suppressed. If one of the detectors 
correctly registers the vacuum, the chance of a single photon being present in the 
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Figure 8.10: Imperfect detectors. Plotted here are (a) the logarithmic negativity Ej^ 
and (b) the entropy 5vN of the states generated after zero (dotted), one (dashed) and two 
(wider dashing) iterations when detectors with efficiency rj are used. The input state is 
\^jJa) defined in the text. 

other detected mode is zero and the chance that a two-photon state is present is only 
1.4%. This means that in a demonstration experiment, if input states generated by a 
Procrustean methods as described in chapter [3 were employed, a large entanglement 
increase could be observed in an implementation of the Gaussification protocol with 
current detector efficiencies. 

While inefficient detectors are clearly a very important experimental imperfection 
for the scheme, there are other obstacles which would make the implementation 
of the full iterative scheme difficult. In the next section we discuss the effect of 
decoherence in the pulses when they are stored between iterations. 



8.4.2 Decoherence between iterations 

For a demonstration experiment of a small number of iterations of the protocol, one 
could avoid the need for storing pulses between iterations by sending through all 
states to their next iterations, and then repeating the experiment many times until 
all measurement outcomes are successful in "one shot" . The probability of success 
will decrease approximately exponentially with increasing number of steps so this 
is not practical for many iterations. To avoid the exponential decrease in success 
probability one would have to store pulses resulting from a successful iteration until 
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Figure 8.11: Imperfect detectors. Here we see how the scheme copes with inefficient 
detection when a mixed input state is used. The degree of entanglement (a) and the 
entropy (b) are plotted against detector efficiency for two iterations as in figure IH^TUl with 
mixed state as input. 
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Figure 8.12: Imperfect detectors. Here we see the trade-off between the mixedness of 
the input state and the necessary efficiency to see an increase in entanglement in the first 
iteration. The logarithmic negativity of the initial state (solid) line and after one iteration 
(dotted line) are plotted as a function of detector efficiency r/ for the state p (defined in 
equations H8.25() to H8.28() ') with A = 1 and r taking the following values: 0.25 (black lines), 
0.5 (blue lines), 0.75 (red lines) and 1.0 (green lines). 

another pair has been successfully generated. Additionally storage would also be 
required for every single pulse generated by a successful iteration, as Alice and Bob 
would need to wait for the signal from each other before they knew if the iteration 
had been successful on both sides. If Alice and Bob are not too far from each other, 
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however, this storage time need not be long, so we shall neglect it here. 

So far, there have been several proposals towards the storage of the quantum 
states of light pulses. The first is to simply divert the light pulse into a loop of 
optical fibre until it is needed. The use of loops of fibre to temporarily store light 
pulses is commonplace in current experiments. However, the ability to switch a 
pulse into and out of storage on demand on the time-scales which would be required 
here is a difficult technological challenge. 

A second, and more radical approach is to use electro-magnetically induced trans- 
parency effects to transfer the quantum state onto collective excitations in atomic 
vapours [FleischhauerOO] . The first experiments on this effect have been promis- 
ing [PhillipsOl], but it remains a complicated procedure, and currently realised pro- 
tocols are suitable only for the storage of qubits, not higher-dimensional states. 
Nevertheless, the storage of quantum states of light pulses is of vital importance 
for many proposals of optical quantum information processing, including quantum 
computation [KnillOl] and is thus a very active research area. 

In this section, therefore, we shall not restrict ourselves to considering a partic- 
ular method of storing light, but shall consider three kinds of generic decoherence 
processes; absorption, dephasing and Gaussian phase defusion. The first of these 
would be particularly important in fibre loop storage and the second could well be 
an important decoherence mechanisms in atomic vapour storage [Mewes04]. The 
third process, phase defusion, would be due to noise in the laser source itself, which 
could mean that pairs of pulses produced at different times are no longer exactly 
phase coherent. 

In order to model these effects, the recurrence relation in equation ()8.13|) needs 
only minor modification for the case where two different states p'-*-' and a'-*-' are used 
as input to the (« + 1) iteration, 

abed 

ya,b,c,d ~ 2^ 2-^ 2-^ 2-^ ^^■^a,b;c,d h' s,t;n,m'-' a-s,b-t;c-n,d-m \O.OOj 
s=0 t=a n=0 m=0 

where the coefficients M^'^:"'J* are given in equation ()8.14|) . Let us assume that one 
of the pair of input states p'-*^ was created directly before the iteration and therefore 
has undergone no decoherence. The other state cr(j) = V[p{i)\ is the result of the 
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Figure 8.13: Absorbing channel. We plot here (a) the logarithmic negativity and (b) the 

entropy of the state p{i) generated by i iterations of the protocol for i = (solid), i = 1 
(dotted), i = 2 (dashed) and i = 3 (wider dashing) where one of the inputs at each stage 
was stored in a Gaussian absorbing channel with transmission coefficient 6. The initial 
state is \ipa)b defined above. 
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Figure 8.14: Absorbing channel. Now we see the effect of the absorbing channel for 
photon storage when the input state in p^- For this input state we plot (a) the entanglement 
and (b) the entropy with teh above conventions. 

same state undergoing decoherence captured by super-operator P. For simplicity we 
assume symmetry between all attempts and also that the same decoherence occurs 
at each iteration. 

The first decoherence process we consider is absorption. 
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8.4.3 Absorbing channel 

Gaussian absorption has already been described in section 17.11 If we set the trans- 
mission coefficient of the storage loops to 6, the channel is described on the matrix 
element level via the following relation 

- E \/(" ] C D ) C I *) ^"''"(^ - 

(8.36) 

Combining equations ()8.35|) and ()8.36p . we simulate two iterations of the protocol 
for the model states {ipa) and \pb)- The results of these simulations for the test-states 
lip) a and pb are plotted in figures IHUHl and IHHH We see that the performance of the 
protocol declines steeply as the transmission coefficient gets lower. As one would 
expect, this is initially accompanied by an increase in the entropy of the states, but 
then the entropy is seen to decline with decreasing r. This is because as r gets 
closer to zero, one of the input states becomes closer to the pure vacuum state. We 
see similar behaviour for the pure input state \tpa) and pb- This is in contrast to the 
case of inefficient detectors where mixed input states were increasingly less robust 
than pure inputs. 

Optical fibres have a high transmission coefficient. For example, infra-red optical 
fibre with an attenuation coefficient of 3dB/km is commercially available [Tec04]. 
This means that the transmission coefficient of a fibre of d km is 0.93"^. For a 
reasonable increase in entanglement, figures 1^. 131 and 18.141 show that the absorbing 
channel with a transmission coefficient of greater than 95% would be desirable. 
This corresponds to a length of fibre of around 750m, or a maximum storage time 
of the order of microseconds. With a pulse repetition rate on the order of GHz and, 
remembering that the Gaussification iterations typically have success probabilities 
on the order of 1/2, one would, if perfect switching were available, be able to store 
pulses for the necessary waiting time between iterations without enough absorption 
occurring to disrupt the protocol. 
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Figure 8.15: Partially Dephasing Channel. Here see the effect of using a partially 
dephasing channel with parameter k on the entanglement (a) and entropy (b) for pure 
input \ipa) ~ Same conventions as figure ISTT^ 

8.4.4 Partially dephasing channel 

Partial dephasing may be an important decoherence process for quantum states of 
light trapped in atomic vapour [Mewes04] or other novel storage methods. It can be 
considered in one of the following equivalent ways; that with probability k, a random 
phase rotation is applied to the state, or that with probability k the environment 
"measures" the photon number of the state. The effect of this is to reduce all off- 
diagonal elements by factor k. Thus it is particularly simple to simulate. When 
both modes of a two-mode state pass through such a channel, the matrix elements 
transform, 

^al,c,d = K + k{1- K)6a,c + " l^)h,d + + ^a,c^b,d (8-37) 

where 5ij = 1 when i = j otherwise 6ij = 0. 

Results of the simulations on input states \psia) and pb are plotted in figures EUHl 
and lOH 



8.4.5 Phase diffusion channel 

In this section we consider the effect of a random decorrelation in the phase of the 
laser pulse. As mentioned briefly in section ITHj noise in the lasing process manifests 
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Figure 8.16: Partially Dephasing Channel. Here see the effect of using a partially 
dephasing channel with parameter k on the entanglement (a) and entropy (b) for mixed 
input ph - Same conventions as figure I5TT^ 
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Figure 8.17: Here we plot the logarithmic negativity for three iterations of the protocol 
on input \ip), where, before each iteration both states pass through a partially dephasing 
channel with dephasing parameter k. 

itself as weak random fluctuations in the phase of the laser source. Thus, the phase 
of the earlier stored pair of pulses, may not be the same as its partner pair in 
the next iteration. We model this as phase diffusing channel acting on the stored 
pulses. We model the phase fluctuations here as a random phase rotation through 
angle with obeying a Gaussian probability distribution with standard deviation 
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Figure 8.18: Phase Diffusion Channel. Here see the effect of storing one pulse in a 
phase diffusing channel with standard deviation v between iterations on the entanglement 
(a) and entropy (b) for pure input \ipa) ~ Same conventions as figure IHTTHl 

V. Furthermore, since both sides of the state were generated at the same time, we 
assume that both suffer an equal phase shift. This channel transforms the matrix 
elements of p*^*-* as follows. 





' 1 02 - 


1 dcj) 


1=C 2^ 






{a + b- 




^ Pa,b,c,d 



i<j){a+b-c-d) (i) 
^ ya,b,c,d 



(8.38) 



The effects of this "phase diffusion" on the Gaussification were, again, simulated 
numerically for input states \ipa) and pb and the results are plotted in figures IHUHl 
and 18.191 Note that, unlike the other two cases, this kind of "decoherence" can be 
prevented by simply making relative phase measurements on the reference pulses 
for each pulse and correcting accordingly. Of course, this would greatly increase the 
complexity of the implementation. 



8.4.6 Comparison of decoherence effects 

As the results of simulations of all three decoherence processes share many similar 
features we shall discuss them together. Firstly, we notice that in all of these cases, 
the decoherence has a much more serious effect on the performance of the protocol 
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Figure 8.19: Phase Diffusion Channel. Here see the effect of storing one pulse in a 
phase diffusing channel with standard deviation v between iterations on the entanglement 
(a) and entropy (b) for mixed input pf, - Same conventions as figure ESI 

that we saw for imperfect detectors. In each step, we see a sharp decrease in the 
entanglement gained for even a small amount of decoherence. We also notice that, 
in contrast to the previous case, the degradation caused to the entanglement gain 
is similar for both input states {ipa) and pb for all three types of decoherence. This 
is an indication that the special symmetry of \ipa) does not shield it from the effect 
of the decoherence processes. 

In the example we have chosen, the input at each stage is asjTiimetric. Since 
the protocol is intended to operate on symmetric pairs of states, we should consider 
whether degradation in the scheme due to these decoherence processes could be more 
due to the asymmetry in the input rather than the decoherence itself. To test this 
we simulated the effect of iterating the procedure on pairs of states which both pass 
through a partially dephasing channel. The results are plotted in figure 18.171 We 
see that the degradation in entanglement for a given dephasing parameter k is much 
worse than the case plotted in figure 18. 151 It is therefore reasonable to conclude 
that it is the decoherence itself, rather than the asymmetry in the input which is 
degrading the protocol in these cases. 

We have presented a variety of decoherence mechanisms which are simple to 
calculate. Of course, in a real situation a combination of these effects and others 
would occur. Nevertheless it is clear that the scheme is not very robust against 
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Figure 8.20: This figure illustrates the problem of mode-matching optical pulses meeting 
on a beam-splitter. If (a) the modes "match" spatially and temporally as well as in 
frequency, then the desired beam splitter operation on the pulses is realised. If (b) the 
modes do not perfectly match, not only is the beam splitter transformation not perfectly 
realised, but the shape of the mode envelope for each pulse is disrupted. 

decoherence between iterations. This emphasises the fact that the development of 
a method of storing the quantum state of light with a low degree of decoherence is 
vital for the experimental implementation of complex optical protocols in quantum 
information science. 

8.4.7 Other potential problems 

The protocol could be implemented in free-space optical system, or in optical fibres. 
While the free space implementation will lead to less absorption, it suffers from 
a technical problem which limits the size of any optical network for laser pulses 
which employs beam splitters. This is the problem of mode matching the pulses in 
incident on each beam splitter, illustrated in figure 111201 Essentially, this means that 
for the desired beam splitter transformation to be obtained between the pulses, they 
must match perfectly, both in frequency and spatial-temporal envelope. In optical 
fibres, fibre couplers replace beam splitters, then the mode-matching is much easier, 
since spatial mode-matching is guaranteed by the fibre. Thus it is expected that 
the full implementation of schemes such as the Gaussification protocol will need a 
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fibre-based implementation. 

8.5 Conclusion 

In this chapter, we have presented an iterative non-deterministic protocol called 
"Gaussification" . The protocol allows the generation of approximately Gaussian bi- 
partite states from a non-Gaussian supply. We have shown that for certain classes 
of input states, the number of iterations need not be too high, before a near Gaus- 
sian form is reached. We have also shown that the entanglement of the states is 
often increased in the protocol. We have comprehensively analysed the properties 
of the scheme when all operations are ideal, and also simulated how the scheme 
may perform when implemented with realistic imperfect operations. We have found 
that the scheme can still perform well when only detectors of currently available 
efficiencies are employed, but that the scheme will be badly degraded if the pulses 
undergo decoherence between iterations. Thus, we believe that while a demonstra- 
tion of entanglement gain in a few iterations of the protocol can be achieved with 
current technology, until an effective means of storing the quantum state of light is 
developed, an implementation with many iterations generating a very highly entan- 
gled state will be out of reach. Nevertheless, even the implementation of a simple 
demonstration experiment would be of enormous importance as it would be the first 
implementation of entanglement distillation in the continuous variable regime and 
a great step forward in our ability to manipulate entanglement in light. 
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Chapter 9 

Summary and Outlook 



The generation of highly entangled states in spatially separated sub-systems promises 
many useful applications, and will allow tests of some of the fundamental proper- 
ties of quantum mechanics. In this thesis, we have presented new proposals for the 
generation and manipulation of entanglement in quantum optical systems. Here we 
present a summary. 

In chapter IH we introduced a proposal for the generation of a maximally en- 
tangled Bell state between two spatially separated cavity modes. The basis of the 
proposal is a very brief interaction between a mediating atom and each cavities. We 
have seen that this gives the proposal a remarkable robustness against experimen- 
tal errors such as imperfect control over the atomic speed and path and detection 
inefficiency. 

In chapter El we described a novel method of creating Bell states in ions trapped 
in spatially separated cavities. This proposal uses continual weak driving of the 
atoms and a detection of a photon after which-path information is erased by a 
beam-splitter. This scheme has many advantages over previous proposals. In the 
strong coupling regime, where decoherence is negligible, and with ideal detectors, 
it becomes near deterministic in spite of the crucial role played by the detectors. 
Furthermore, when likely experimental imperfections are taken into account the 
scheme still allows high fidelity Bell states to be generated. 

In chapter [3 we showed how simple Procrustean linear optical protocols can 
generate highly entangled non-Gaussian states from a supply of weakly entangled 
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Gaussian states, and we introduced, in chapter |Hl an iterative procedure which 
we call "Gaussification" . This is a simple linear optical protocol, which generates, 
from a large supply of non-Gaussian states, a small number of Gaussian states 
which can be significantly more entangled than the input supply. Combining these 
two procedures gives a linear optical entanglement distillation procedure. We have 
assessed the performance of the protocol when important experimental imperfections 
are taken into account. We showed that the protocol can still produce significant 
gains in entanglement when detectors with the efficiencies currently available are 
employed, and thus would be suitable for a demonstration experiment with current 
detectors. We also showed that the scheme can be badly degraded by decoherence 
in the pulses when stored between iterations. This highlights the experimental 
challenge a full many-iteration implementation of the scheme would pose. 

There are many avenues for future work based upon these results. The schemes 
presented in chaptersEJandElused an approach of weak interaction plus measurement 
to allow robust entangled state preparation. We have focussed upon cavity QED 
implementations, but this principle would naturally lend itself to other kinds of 
system. In particular, in systems where couplings exist, efficient measurement may 
be made and robustness is needed, such as coupled charge qubits [PashkinOS], this 
approach could be very fruitful. 

In chapter the detection of a photon induces a projection onto a maximally 
entangled state of the two atoms. It has recently been shown [Browne04] that this 
kind of projection can be used to create special multi-particle entangled states called 
"cluster states" [BriegelOl] which have many interesting properties. In particular, 
they are a resource for quantum computation through local measurements alone 
[RaussendorfOS] . It would be a simple matter to apply these ideas to the setup 
we have outlined, which would then lead to a robust method for generating useful 
multi-party entangled states in the cavity QED setting. 

The entanglement distillation scheme we have presented in chapter 18. 4t while 
promising in its current form, could still be improved. In particular the Procrustean 
entanglement concentration step, which incorporates the non-Gaussian operation for 
the distillation, has a very low success probability. Recent results have suggested 
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that the cross-Kerr non-hnearities one can produce using electro-magnetically in- 
duced transparency in atomic vapours [Schmidt96] , can be useful for the manipula- 
tion of quantum states of light beyond the linear optical regime [Munro03,Barrett04] . 
Such operations are certainly non-Gaussian and it would be worthwhile investigating 
whether a more effective non-Gaussian preparation step might be developed using 
these techniques. 

We have shown that the Gaussification process generates pure entangled states 
from certain particular mixed input states. In principle, all input states could be 
transformed into this special form by an appropriate completely positive map. How- 
ever, this does not mean that either the map is easy to implement experimentally, 
or that the state produced after the resultant Gaussification would still have more 
entanglement than the input. Work on characterising the general operations which 
can be achieved linear optically is progressing [Ahnert04] and it would be very in- 
teresting to apply these ideas to this case. 

The most important obstacle to the full implementation of the Gaussification 
protocol and linear optical quantum computation is that no effective method of 
storing quantum states of light between non-deterministic steps exists. Whilst im- 
portant first steps towards this are being made in atomic vapour experiments [Fleis- 
chhauerOO, PhiUipsOl], this work emphasises the importance of investigating other 
novel kinds of quantum memory for light. 

The measurement of quantum systems has long been our best connection with 
the strange phenomena which occur at quantum scales, and it has played an in- 
creasingly important and diverse role as the field of quantum information science 
has developed. In this thesis we have demonstrated the power that measurement 
can wield in particular for the preparation of useful quantum states and the manip- 
ulation of quantum information. We hope that these ideas can make a contribution 
towards the realisation of quantum information science's extraordinary potential. 
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